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Abstract: We develop several methods that allow us to compute all- loop partition 
functions in perturbative Chern-Simons theory with complex gauge group Gc, sometimes 
in multiple ways. In the background of a non-abelian irreducible flat connection, per- 
turbative Gc invariants turn out to be interesting topological invariants, which are very 
different from finite type (Vassiliev) invariants obtained in a theory with compact gauge 
group G. We explore various aspects of these invariants and present an example where 
we compute them explicitly to high loop order. We also introduce a notion of "arithmetic 
TQFT" and conjecture (with supporting numerical evidence) that SL{2, C) Chern-Simons 
theory is an example of such a theory. 
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4.2.2 



1. Introduction and summary 

Three-dimensional Chern-Simons gauge theory is a prominent example of a topological 
quantum field theory (TQFT). By now, Chern-Simons theory with a compact gauge group 



G is a mature subject with a history going back to the 1980's (see e.g. |16, for excellent 
reviews) and with a wide range of applications, ranging from invariants of knots and 3- 
manifolds [^] on one hand, to condensed matter physics |25, ^] and string theory |^5| on 
the other. 

In this paper, we will be interested in a version of Chern-Simons gauge theory with 
complex gauge group Gc. Although at first it may appear merely as a variation on the 
subject, the physics of this theory is qualitatively different from that of Chern-Simons 
gauge theory with compact gauge group. For example, one important difference is that to 
a compact Riemann surface S Chern-Simons theory with compact gauge group associates 
a finite-dimensional Hilbert space TCj:, whereas in a theory with non-compact (and, in 
particular, complex) gauge group the Hilbert space is infinite-dimensional. Due to this and 
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other important differences that wih be explained in further detail below, Chern-Simons 
gauge theory with complex gauge group remains a rather mysterious subject. First steps 



toward understanding this theory were made in [53| and, more recently, in |22]. 

As in a theory with a compact gauge group, the classical action of Chern-Simons gauge 
theory with complex gauge group Gc is purely topological — that is, independent of the 
metric on the underlying 3-manifold M. However, since in the latter case the gauge field A 
(a gc-valued 1-form on M) is complex, in the action one can write two topological terms, 
involving A and A: 

S = — [ Tr(AAdA+-AAAAA) (1.1) 
8vr Jm V 3 / 

+ — / Trf^ A + A ^ A ^ 
Svr Jm V 3 

Although in general the complex coefficients ("coupling constants") t and t need not be 
complex conjugate to each other, they are not entirely arbitrary. Thus, if we write t = k + a 
and i = k—a, then consistency of the quantum theory requires the "level" k to be an integer, 
A; S Z, whereas unitarity requires a to be either real, o" S M, or purely imaginary, a S iM; 
see e.g. [||]. 

Given a 3-manifold M (possibly with boundary), Chern-Simons theory associates to 
M a "quantum Gc invariant" that we denote as Z{M). Physically, Z[M) is the partition 
function of the Chern-Simons gauge theory on M, defined as a Feynman path integral 

Z{M) = I e}^ VA (1.2) 



with the classical action (1^1). Since the action ( |l.lj ) is independent of the choice of metric 
on M, one might expect that the quantum Gc invariant Z{M) is a topological invariant 
of M. This is essentially correct even though independence of metric is less obvious in the 
quantum theory, and Z[M) turns out to be an interesting invariant. How then does one 
compute Z[M)1 

One approach is to use topological invariance of the theory. In Chern-Simons theory 
with compact gauge group G, the partition function Z{M) can be efficiently computed by 
cutting M into simple "pieces," on which the path integral ( p^ ) is easy to evaluate. Then, 
via "gluing rules," the answers for individual pieces are assembled together to produce 
Z[M). In practice, there may exist many different ways to decompose M into basic building 
blocks, resulting in different ways of computing Z{M). 

Although a similar set of gluing rules should exist in a theory with complex gauge 
group Gc, they are expected to be more involved than in the compact case. The underlying 
reason for this was already mentioned: in Chern-Simons theory with complex gauge group 
the Hilbert space is infinite dimensional (as opposed to a finite-dimensional Hilbert space 
in the case of compact gauge group G). One consequence of this fact is that finite sums 
which appear in gluing rules for Chern-Simons theory with compact group G turn into 
integrals over continuous parameters in a theory with non-compact gauge group. This is 
one of the difficulties one needs to face in computing Z{M) non-perturbatively, i.e. as a 
closed-form function of complex parameters t and i. 
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A somewhat more modest goal is to compute Z{M) perturbatively, by expanding the 
integral (^]^) in inverse powers of t and i around a saddle point (a classical solution). In 
Chern-Simons theory, classical solutions are flat gauge connections, that is gauge connec- 
tions A which obey 

dA + A^A = 0, (1.3) 

and similarly for A. A flat connection on M is determined by its holonomies, that is by a 
homomorphism 

p:^i(M)^Gc. (1.4) 

Of course, this homomorphism is only defined modulo gauge transformations, which act 
via conjugation by elements in Gc- 

Given a gauge equivalence class of the flat connection A^ or, equivalently, a conjugacy 
class of the homomorphism /j, one can define a "perturbative partition function" Z^p") (M) 



by expanding the integral (1.2) in inverse powers of t and i. Since the classical action (|L 
is a sum of two terms, the perturbation theory for the fields A and A is independent. As a 
result, to all orders in perturbation theory, the partition function Z^f\M) factorizes into 
a product of "holomorphic" and "antiholomorphic" terms: 

Z^p\M) = Z^P\M;t)Z'^P\M;t). (1.5) 

This holomorphic factorization is only a property of the perturbative partition function. 
The exact, non-perturbative partition function Z(M) depends in a non-trivial way on both 
t and t, and the best one can hope for is that it can be written in the form 

Z{M) = ^Z^P\M;t)Z^P\M;t), (1.6) 



where the sum is over classical solutions ( |1.3D or, equivalently, conjugacy classes of homo- 
morphisms ( |1.4D . 

In the present paper, we study the perturbative partition function Z^p\M). Due to the 
factorization (1^), it suffices to consider only the holomorphic part Z^P\M;t). Moreover, 



since the perturbative expansion is in the inverse powers of i, it is convenient to introduce 
a new expansion parameter 

/i = — , (1.7) 

which plays the role of Planck's constant. Indeed, the semiclassical limit corresponds to 
— > 0. In general, the perturbative partition function Z^p\M; K) is an asymptotic power 
series in h. To find its general form one applies the stationary phase approximation to the 
integral ([L^): 



Z^P\M;h) 



\ n=0 / 



This is the general form of the perturbative partition function in Chern-Simons gauge 
theory with any gauge group, compact or otherwise. It follows simply by applying the 
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stationary phase approximation to the integral (|1.2| ), which basicahy gives the standard 
rules of perturbative gauge theory |52, ^, ^, S] that will be discussed in more detail below. 



For now, we note that the leading term is the value of the classical Chern-Simons 
functional evaluated on a flat gauge connection A^^^ associated with a homomorphism p, 

= -7 / Trf^^") A d^^") + Ia^p'^ a A A'^pA . (1.9) 
4 Jm V 3 / 

Moreover, the coefficient of the next-leading term, 5'^p\ is an integer (which, like all other 
terms, depends on the 3-manifold M, the gauge group Gc, and the classical solution p). 
Each coefficient Sn"^ is obtained by summing over Feynman diagrams with n loops. For 
example, the "two-loop term" 5*2^^ is obtained by adding contributions of the two kinds of 
Feynman diagrams shown in Figure |l[ Although certain features of the Feynman rules in 
Chern-Simons theory with complex gauge group Gq are similar to those in a theory with 
compact gauge group G, we will see that there exist important differences. 





a) b) 

Figure 1: Two kinds of 2-loop Feynman diagrams that contribute to S'2'''' in Chern-Simons gauge 
theory. 

In Chern-Simons theory with compact gauge group, one often develops perturbation 
theory in the background of a trivial (or reducible) flat connection A^^\ As a result, the 
perturbative coefficients Sn^ have a fairly simple structure; they factorize into a product 
of topological invariants of M — the finite type (Vassiliev) invariants and varations thereof 



and group theory factors [33|. In particular, they are rational numbers. 



The situation is very different in Chern-Simons theory with complex gauge group; the 
are completely new, unexplored invariants of M. From the physics point of view, one 
novelty of the complex case is that it involves perturbation theory in the background of a 
genuinely non-abelian flat connection, a problem that has not been properly addressed in 
the existing literature. One of our goals below will be to make a few modest steps in this 
direction. As we shall see, in general, the information about the gauge group Gc and the 
3-manifold M is mixed within (M) in a non-trivial way. In particular, in this case the 
s!^\Mys are not finite type invariants and are typically not valued in Q. 

If we take p to be the trivial representation, denoted as p = (and corresponding 
to a trivial flat connection A'^^'^ = 0), the perturbative invariants Sn\M) are actually 
identical to those in Chern-Simons theory with compact gauge group G. Then Z^^\M;h) 
can be obtained from the perturbative partition function of Chern-Simons theory with 
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gauge group G simply by allowing h to be complex,^ 



zP(M-n) = Z§\M;h). (1.10) 



Both sides of this equation have the form of the asymptotic series ( |1.8| ) with S^^ = 0. 

This example, however, is very special, as one can see e.g. by considering hyperbolic 
3-manifolds, which, in a sense, constitute the richest and the most interesting class of 3- 
manifolds. A hyperbolic structure on a 3-manifold M corresponds to a discrete faithful 
representation of the fundamental group 7ri(M) into Isom^(]H['^) = PSL{2, C), the group of 
orientation-preserving isometrics of 3-dimensional hyperbolic space M^. One can compose 
this representation with a morphism (p from PSL(2) to a larger algebraic group to obtain 
a group homomorphism 

p:TTi{M)^ PSL{2,C) ^Gc. (1.11) 

The flat connection A^''^ associated with this homomorphism is genuinely non-abelian and 
the corresponding perturbative invariants are interesting new invariants of the hyper- 
bolic 3-manifold M. See Table § on page ^ for the simplest example of this type. 

Standard textbooks on quantum field theory tell us how to systematically compute 
the perturbative invariants S'if'* for any n, starting with the classical action (^^) and a 
classical solution (^]^), and using the Feynman diagrams mentioned above. This compu- 
tation will be reviewed in section 2 for the problem at hand. However, direct computation 
of Feynman diagrams becomes prohibitively complicated when the number of loops n be- 
comes large. Therefore, while Feynman diagrams provide the first-principles approach to 
calculating Z^p\M;h), it is useful to look for alternative ways of computing and defining 
the perturbative Gc invariants of M. 

Such alternative methods often arise from equivalent physical descriptions of a given 
system. In some cases, interesting equivalent descriptions are obtained merely via differ- 
ent gauge choices. A relevant example is the much-studied Chern-Simons gauge theory 
with compact gauge group G. Perturbative expansion of this theory around the trivial 
flat connection A^^^ = in the covariant Landau gauge leads to the configuration space 
integrals for the Vassiliev invariants, while in the non-covariant light-cone gauge it leads 
to the Kontsevich integral [^]. This example illustrates how a simple equivalence of two 
physical descriptions (in this case, based on different gauge choice) leads to an interesting 
mathematical statement. Much deeper examples involve non-trivial duality symmetries in 
quantum field theories. For example, duality was the key element in the famous work of 



Seiberg and Witten [45| that led to an exact description of the moduli space of vacua in 
M = 2 gauge theory and to a new way of computing Donaldson invariants of 4-manifolds. 

Similarly, in the present case of Chern-Simons gauge theory with complex gauge group, 
there exist several ways to quantize the theory. When combined with topological invari- 
ance, they lead to multiple methods for computing the perturbative invariant Z^''^(M; h). 
Mathematically, each of these methods can be taken as a definition of Z^p\M; h); however. 



^In a theory with compact gauge group, one frequently uses the "level" k as the coupling constant. It is 
related to h via ft — iTi/{k + ft^), where is the dual Coxeter number of G, e.g. = N for G = SU{N). 
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it is important to realize that they all have a common origin and are completely equivalent. 
In the mathematical literature, this is usually stated in the form of a conjecture. To wit, 
the following are equivalent: 



Feynman diagrams: This is the most traditional approach, which was already 



mentioned earlier and which will be discussed in more detail in Section 2.1. It gives 
as an asymptotic power series of the form ( |1.8p where every coefficient 
Sn can be systematically computed by evaluating Feynman diagrams. 

Quantization of 7Wflat(Gc, S): According to the general axioms of topological quan- 
tum field theory 0, Chern-Simons gauge theory associates a Hilbert space TIy, to a 
closed surface S (not necessarily connected), and, similarly, it associates a vector in 
to every closed 3-manifold M with boundary S. In Chern-Simons theory with 
complex gauge group Gc, the space TLy. is obtained by quantizing the classical phase 
space A4flat(Gc, S), which is the moduli space of flat Gc connections on S. Therefore, 
in this approach Z^P'f (M; K) is obtained as the wave function of a state associated with 
M. (If M is a closed 3-manifold without boundary, one can consider e.g. a Heegaard 
decomposition M = M_|_ Us along S and apply this approach to each handlebody, 
M+ and M_.) In particular, as we explain in Section |2.2| , Z^p\M; K) obeys a system 
of Schrodinger-like equations 

Ai Z^P\M-h) = 0, (1.12) 

which, together with appropriate boundary conditions, uniquely determine Z^p\M; h). 
One advantage of Chern-Simons theory with complex gauge group is that the classical 
phase space 7Wfiat(G'c, S) is a hyper-Kahler manifold, a fact that considerably sim- 
plifies the quantization problem in any of the existing frameworks (such as geometric 
quantization [^] , deformation quantization |6|, |3^ , or "brane quantization" [^] ) . 

"Analytic continuation": Since Gc is a complexification of G, one might expect 
that Chern-Simons theories with gauge groups G and Gc are closely related. This is 
indeed the case |22]. For example, if M = S'^ \ X is the complement of a knot (or 
a link) in the 3-sphere, one can compare the perturbative Gc invariant Z^p^ {M; h) 
with the expectation values of a Wilson loop supported on C in Chern-Simons 
theory with compact gauge group G. The definition of the Wilson loop operator, 
Wr{K) := Tvji Pex.p<fj^A, involves a choice of a representation R of G, and the 
expectation value 

PG,R{K;q) := {Wr{K)) (1.13) 

is a certain polynomial invariant of K. Letting R = Rx be the irreducible representa- 
tion of G with highest weight A, the invariant quadratic form — Tr on the Lie algebra 
Q = Lie(G) identifies A with an element A* of the Cartan subalgebra t C g. Then, 
one considers a double-scaling limit such that 

h^O, X*^oo, u := hX* = hxed . (1.14) 
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Taking this limit and analytically continuing to complex values of u, one has [22| 
Z^p\M;h) 



asymptotic expansion of Pg,r^{K; q) , (1-15) 



where q = e and p is determined by u. Further details will be explained in Sec- 



tion 2.2.1. 



State sum model: Finally, using the topological nature of the theory, one can 
construct Z^^^M; h) by decomposing M into more elementary pieces. In particular, 
one can consider a decomposition of M into tetrahedra, i.e. a triangulation. 



N 



M = IjAi. (1.16) 



i=l 



Since each tetrahedron Aj is a 3-manifold with boundary, it defines a state (a vector) 
in the Hilbert space associated to the boundary of Aj. Then, gluing the tetrahedra 
together as in (P^ ) leads to an expression for Z(M) as a sum over different states. 



Of course, the result cannot depend on the choice of the triangulation. A well-known 



example of such a state sum model is the Turaev-Viro invariant [50|. 

In general, states associated to the tetrahedra Aj are labeled by representations of the 
gauge group, which means that the set of labels is discrete in a theory with compact 
gauge group, but is continuous in a theory with non-compact gauge group. Therefore, 
in Chern-Simons theory with complex gauge group Gc one might hope to define the 
perturbative Gc invariant Z^P\M]h) as an integral over states, starting from the 
initial data of a triangulation ( |1.16| ) and a homomorphism p. We will therefore use 
the word "state integral" instead of using the more familiar term "state sum" for 
invariants defined in this way. Constructing such an invariant for the simplest theory 
with Gc = SL(2,C) is one of the main goals of the present paper, and will be 
considered in detail in section ^ 

Of the different approaches just outlined, the first three have been used previously to tackle 
Chern-Simons theory with complex gauge group, while the fourth is new. In principle one 
can compute Z^p'^ (M; K) following any of these four approaches individually. However, 
as we show in Section |2|, combining them together dramatically improves computational 
power. For example, when M is a knot complement, one can combine approaches 2 and 
3 to learn that the operators Ai in ( |1.12 ) should also annihilate the polynomial invariants 



PG,R{K;q) computed by Chern-Simons gauge theory with compact gauge group G, 

Ai PG,R{K;q) = 0. (1.17) 
(The normalization by ^(S^) is not important here.) This provides a very simple method 



for finding the Aj's — see sections and 4.2.1 for further details and a concrete example. 



In turn, combining approaches 2 and 3 with the first approach implies that, in the 
limit ( |1.14| ), the asymptotic expansion of the polynomial invariant Pg,r{K; q) has the form 
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For G = SU{2) (whose complexification is Gq = SL{2,C)) this statement is aheady 
known as a generahzation of the Volume Conjecture, first proposed in [^] and further 
refined in In particular, if G = SU{2) and R is its n-dimensional representation, 

the polynomial Pg,r{K; q) is called the "n-colored Jones polynomial" of K and is denoted 
Jn{K;q). The proposed relation ( 1.15| ) then says that the asymptotic expansion of the 



colored Jones polynomial reproduces the perturbative 5L(2,C) invariant Z^P\M;h); in 
particular, the expansion has the structure (1.8) with a very interesting dependence on 
both p and h. 



In Section we describe the fourth approach. Building on the work of Hikami |26 
27], we construct a state sum model for perturbative SL{2,C) Chern-Simons theory on a 
hyperbolic 3-manifold M (possibly with boundary). Schematically, the state sum model 
has the form 

„ N N-bo{T) 

JCp viirn 

where ^h{Aj) is the quantum dilogarithm function, associated to every tetrahedron Aj in 
a triangulation of M, 6o($]) is the number of connected components of S, and /(p) is a 
simple quadratic polynomial of the integration variables pi. 

One limitation of the construction of the state sum model, in the form stated here, 
is that it applies only to hyperbolic 3-manifolds. Although we expect that it can be 
extended to arbitrary 3-manifolds with boundary, we will not presently attempt to seek 
such a generalization. Hyperbolic 3-manifolds do provide the richest and, in a sense, the 
most interesting laboratory for studying SL{2, C) Chern-Simons theory. There are various 
other questions that we do not address here. One obvious example is to define a state sum 
model for perturbative Chern-Simons theory with complex gauge group Gc not equal to 
SL{2,C) (or PSL{2,C) = 5'L(2, C)/{±1}). A much deeper issue would be to go beyond 
perturbation theory, even in the simplest case of Gc = SL{2,C). Indeed, once we have a 
state sum model that allows the computation of all perturbative invariants s!^^ to arbitrary 
order n, the next natural question is whether there exists a state sum model that allows a 
non-perturbative computation of Z{M), as a function of q = e^^. We leave these questions 
for future work. 

By the end of this paper we will have developed several tools that allow us to solve 
perturbative Chern-Simons gauge theory with complex gauge group exactly, in some cases 
in more than one way. By an exact solution we mean a straightforward method that 
explicitly determines (at least in principle) all perturbative invariants Sn^ . In Section |^, 
we illustrate how these methods work in the case of the figure-8 knot complement (the 
simplest hyperbolic knot complement) and compare the results of different calculations. 
As expected, we find perfect agreement. 

Even in this simplest example, the perturbative invariants turn out to be highly 
non-trivial (see e.g. Table |2|). As we illustrate in this example, the new perturbative 
invariants Sn^ computed by Chern-Simons gauge theory with complex gauge group Gc in 
the background of a non-abelian irreducible flat connection A^''^ contain much interesting 
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information about the geometry of M. We believe that studying these invariants further, 
both physically and mathematicahy, can lead to exciting developments in both subjects. 



2. Perturbation theory around a non-trivial flat connection 

In this section we pursue traditional approaches to Chern-Simons gauge theory based on the 
evaluation of Feynman diagrams and the quantization of moduli spaces of flat connections. 
In particular, we explain how one can systematically compute the perturbative Gc invariant 
Z'^P^M; h) to arbitrary order in the ^-expansion for non-trivial flat connections. Through 
the analysis of Feynman diagrams, we will find that the coefficients Sn^ have a very special 
structure, motivating us to define the notion of "arithmetic QFT." 

2.1 Feynman diagrams and arithmetic TQFT 



Let us examine more carefully each term in the perturbative expansion (1.8) 

= exp 




1 °° \ 

-5(^)log;.+ ^5S,n" . (2.1) 

n=0 / 



As already mentioned in the introduction, the leading term is the value of the classical 
Chern-Simons functional ( |1.9D evaluated on a flat gauge connection A^^'^ that corresponds 
to a homomorphism p : 7ri(M) — > Gc- It is also easy to understand the integer coefficient 
in (|2.lD . The homomorphism p defines a flat Gc bundle over M, which we denote as 
Ep. Letting H^{M; Ep) be the i-th cohomology group of M with coefficients in the flat 
bundle Ep, the coefficient 6^^^ is given by 

= h^-h\ (2.2) 

where /i* := dim W (M ; Ep). Both this term and s[^^ come from the "one- loop" contribu- 
tion to the path integral ([L.2|); s[^^ can be expressed in terms of the Ray-Singer torsion of 
M with respect to the flat bundle Ep {cf. I, g ) , 

S<.).^,„,(^). p.3) 

The geometric interpretation of the higher-order terms Sn^ with n > 1 is more interest- 
ing, yet less obvious. To understand it better, we note that the saddle point approximation 
to the path integral (|l.2|) gives an expression for Sn ^ as a sum of Feynman diagrams with 



n 



loops. Since the Chern-Simons action (1.1) is cubic, all vacuum diagrams (that is, Feyn- 



man diagrams with no external lines) are closed trivalent graphs; lines in such graphs have 
no open end-points. Therefore, a Feynman diagram with n loops (n > 1) has 3(n — 1) 
line segments with end-points meeting at 2(n — 1) trivalent vertices. Each such diagram 
contributes an integral of the form (c/. P]) 



L" , (2.4) 

JV/2n-2 
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where M^"~^ denotes a product of 2n — 2 copies of M and denotes a wedge product 

of a 2-form L{x,y) G Q?'{Mx x Mj^;0c ^SSc)- The 2-form L{x,y), called the "propagator," 
is a solution to the first-order PDE, 

dj,(p)L{x,y) = d^{x,y), (2.5) 

where 5^{x,y) is a J-function 3-form supported on the diagonal in x My and dj^(p) is 
the exterior derivative twisted by a flat connection A^^^^ on the Gc bundle Ep. 

For example, suppose that M is a geodesically complete hyperbolic 3-manifold of finite 
volume. As we pointed out in the introduction, such 3-manifolds provide some of the most 
interesting examples for Chern-Simons theory with complex gauge group. Every such M 
can be represented as a quotient 

M = n^/V (2.6) 

of the hyperbolic space by a discrete, torsion-free subgroup F C P5L(2,C), which is a 
holonomy representation of the fundamental group 7ri(M) into Isom"*'(EI^) = PSL{2,C). 
In what follows, we call this representation "geometric" and denote the corresponding flat 
P5L(2,C) connection as ^fe'^""^). 

In order to explicitly describe the flat connection A.^s^"^\ recall that can be defined 
as the upper half-space with the standard hyperbolic metric 



1 



ds^ = —{dxi + dxi + dxi), xs > . (2.7) 

The components of the vielbein and spin connection corresponding to this metric can be 
written as 

I dxi I dx2 

X3 X3 



dx2 2 dxi 

X3 X3 



3 dx3 3 

X3 

These satisfy de" + e'''"'wt A ec = 0. It is easy to check that the corresponding PSL{2, C) 
connection 

^fe--) = -{w + ^e) = -1- f'^^^ ^dx^-2^dx2\ 

2x3 \ -dx3 



is flat, i.e. obeys (|1.3|). In Chern-Simons theory with gauge group Gq = PSL{2,C), this 
gives an explicit expression for the flat gauge connection that corresponds to the hyperbolic 
structure on M. In a theory with a larger gauge group, one can also define a geometric 
connection ^(s"^"™) by embedding ( |2.8| ) into a larger matrix with extra rows and columns 
of zero entries. We note, however, that ^(8*^°™) constructed in this way is not unique and 
depends on the choice of embedding, cf. (|1.11 ). Nevertheless, we will continue to use the 



notation even in the higher-rank case whenever the choice of embedding introduces 

no confusion. 
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The action of Ton can be conveniently expressed by identifying a point (xi, X2, 2:3) G 
with a quaternion q = xi + X2i + X2,j and defining 



^:g^(ag + 6)/(cg + d), 7=1^" ^jGP5L(2,C). (2.9) 

Exphcitly, setting z = xi + 2x2, we find 

7(^ + 2;3i) = ^' + 4i > (2-10) 

where 



/ ^ {az + b){cz + d) + acxl ^, ^ X3 

~ \cz + d\'^ + \c\'^xl ' ^ Icz + dp + Icpxg 

Let Lo{x,y) be the propagator for M^, i.e. a solution to equation (2^) on x 
with the non-trivial flat connection ^(s°°™). Then, for a hyperbolic quotient space (p.6|), 
the propagator L(x, y) can simply be obtained by summing over images: 

L{x,y) = ^Lo{x,^y). (2.12) 

7er 

Now, we would like to consider what kind of values the perturbative invariants 
can take. For n > 1, the sii'^'s are given by sums over Feynman diagrams, each of which 
contributes an integral of the form (^|^). A priori the value of every such integral can 
be an arbitrary complex number (complex because we are studying Chern-Simons theory 
with complex gauge group) that depends on the 3-manifold M, the gauge group Gq, and 
the classical solution p. However, for a hyperbolic 3-manifold M = H^/F and for the 
flat connection A^seom) associated with the hyperbolic structure on M, we find that the 
^^geom),^ are significantly more restricted. 

Most basically, one might expect that the values of 5i^*^°™'''s are periods 

^feeom) g p (2.13) 

Here V is the set of all periods, satisfying 

QcQcPcC. (2.14) 

By definition, a period is a complex number whose real and imaginary parts are (absolutely 
convergent) integrals of rational functions with rational coefficients, over domains in M" 
defined by polynomial inequalities with rational coefficients |31]. Examples of periods are 



powers of vr, special values of L-functions, and logarithmic Mahler measures of polynomials 
with integer coefficients. Thus, periods can be transcendental numbers, but they form a 
countable set, V. Moreover, V is an algebra; a sum or a product of two periods is also a 
period. 

Although the formulation of the perturbative invariants S^^'^^ in terms of Feynman 



diagrams naturally leads to integrals of the form (|2.4| ), which have the set of periods, V, 
as their natural home, here we make a stronger claim and conjecture that for n > 1 the 
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^^geom),^ are algebraic numbers, i.e. they take values in Q. As we indicated in (2.14), the 
field Q is contained in V, but leads to a much stronger condition on the arithmetic nature 
of the perturbative invariants sl^^°^\ In order to formulate a more precise conjecture, we 
introduce the following definition: 



Definition: A perturbative quantum field theory is called arithmetic if, for all n > 1, the 
perturbative coefficients take values in some algebraic number field K, 

G K, (2.15) 

and S[''^ G Q-logK. 

Therefore, to a manifold M and a classical solution p an arithmetic topological quan- 
tum field theory (arithmetic TQFT for short) associates an algebraic number field K, 

{M,p)-^K. (2.16) 

This is very reminiscent of arithmetic topology, a program proposed in the sixties by 
D. Mumford, B. Mazur, and Yu. Manin, based on striking analogies between number 
theory and low-dimensional topology. For instance, in arithmetic topology, 3-manifolds 
correspond to algebraic number fields and knots correspond to primes. 

Usually, in perturbative quantum field theory the normalization of the expansion pa- 
rameter is a matter of choice. Thus, in the notations of the present paper, a rescaling 
of the coupling constant h ^ Xh hy a numerical factor A is equivalent to a redefnition 
S'n''^ X^~'^S^\ While this transformation does not affect the physics of the perturbative 
expansion, it certainly is important for the arithmetic aspects discussed here. In particular, 
the above definition of arithmetic QFT is preserved by such a transformation only if A G Q. 
In a theory with no canonical scale of h, it is natural to choose it in such a way that makes 
the arithmetic nature of the perturbative coefficients as simple as possible. However, 
in some cases (which include Chern-Simons gauge theory), the coupling constant must obey 
certain quantization conditions, which, therefore, can lead to a "preferred" normalization 
of h up to irrelevant Q- valued factors. 

We emphasize that our definition of arithmetic QFT is perturbative. In particular, it 
depends on the choice of the classical solution p. In the present context of Chern-Simons 
gauge theory with complex gauge group Gc, there is a natural choice of p when M is 
a geodesically complete hyperbolic 3-manifold, namely the geometric representation that 
corresponds to ^(8"^°™). In this case, we conjecture: 



Conjecture 1 (Arithmeticity conjecture): As in ( |2.6D , let M be a geodesically com- 
plete hyperbolic 3-manifold of finite volume, and let p = geom be the corresponding discrete 
faithful representation of tii{M) into PSL{2,C). Then the perturbative Chern-Simons the- 
ory with complex gauge group Gq = PSL{2,C) (or its double cover, SL{2,C)) in the 
background of a non-trivial flat connection A^s^om) arithmetic on M . 
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In fact, we can be a little bit more specific. In all the examples that we studied, we 
find that, for M as in e and for all values of n > 1, the perturbative invariants 5^^°^) 



take values in the trace field of F, 

^(geom) g Q(trr), (2.17) 

where, by definition, Q(tr F) is the minimal extension of Q containing tr 7 for all 7 G F. We 
conjecture that this is the case in general, namely that the SL{2, C) Chern-Simons theory 
on a hyperbolic 3-manifold M = H^/F is arithmetic with K = Q(trF). This should be 
contrasted with the case of a compact gauge group, where one usually develops perturbation 
theory in the background of a trivial flat connection, and the perturbative invariants Sn^ 
turn out to be rational numbers. 

Of course, in this conjecture it is important that the representation p is fixed, e.g. by 
the hyperbolic geometry of M as in the case at hand. As we shall see below, in many cases 
the representation p admits continuous deformations or, put differently, comes in a family. 
For geometric representations, this does not contradict the famous rigidity of hyperbolic 
structures because the deformations correspond to incomplete hyperbolic structures on M. 
In a sense, the second part of this section is devoted to studying such deformations. As 
we shall see, the perturbative Gc invariant Z^P\M;h) is a function of the deformation 
parameters, which on the geometric hrancf? can be interpreted as shape parameters of the 
associated hyperbolic structure. 

In general, one might expect the perturbative coefficients S^^^ to be rational functions 
of these shape parameters. Note that, if true, this statement would imply the above 
conjecture, since at the point corresponding to the complete hyperbolic structure the shape 
parameters take values in the trace field Q(trF). This indeed appears to be the case, at 
least for several simple examples of hyperbolic 3-manifolds that we have studied. 

For the arithmeticity conjecture to hold, it is important that h is defined (up to Q- 



valued factors) as in (1.7), so that the leading term 5g''^ is a rational multiple of the classical 



Chern-Simons functional, c/. (|1.9| ). This normalization is natural for a number of reasons. 
For example, it makes the arithmetic nature of the perturbative coefficients Sn^ as clear 
as possible. Namely, according to the above arithmeticity conjecture, in this normalization 
^(gcom) ^ period, whereas 5*,^^°™^ take values in Q. However, although we are not going 
to use it here, we note that another natural normalization of h could be obtained by a 
redefinition h ^ Xh with A G (27rz)~^ • Q. As we shall see below, this normalization 
is especially natural from the viewpoint of the analytic continuation approach. In this 
normalization, the arithmeticity conjecture says that all S^^°^^ are expected to be periods. 
More specifically, it says that Sn^ G (27ri)"~^ • Q, suggesting that the n-loop perturbative 
invariants s!/^^°^^ are periods of (framed) mixed Tate motives Q(n — 1). In this form, the 
arithmeticity of perturbative Chern-Simons invariants discussed here is very similar to the 
motivic interpretation of Feynman integrals in ||7| . 

Finally, we note that, for some applications, it may be convenient to normalize the 
path integral (|1.2|) by dividing the right-hand side by Z{S^;h). (Since 7ri(S^) is trivial. 



^i.e. on the branch containing the discrete faithful representation. 
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we have Z{S'^; h) = Z^^\S^; h).) This normahzation does not affect the arithmetic nature 
of the perturbative coefficients S'i^'^"™^ because, for M = S^, all the Sn'^'s are rational 
numbers. Specifically, 

^(«'^*) = (^)'"(^)"'n^-'>w»-.)). (2.18) 

aeA+ 

where the product is over positive roots a G A^, r is the rank of the gauge group, and g is 
half the sum of the positive roots, familiar from the Weyl character formula. Therefore, in 



the above conjecture and in eq. ( 2.17 ) we can use the perturbative invariants of M with 
either normalization. 

The arithmeticity conjecture discussed here is a part of a richer structure: the quan- 
tum Gc invariants are only the special case x = of a collection of functions indexed by 
rational numbers x which each have asymptotic expansions in h satisfying the arithmetic- 
ity conjecture and which have a certain kind of modularity behavior under the action of 
S'L(2, Z) on Q |59|. A better understanding of this phenomenon and its interpretation will 
appear elsewhere. 

2.2 Quantization of A4flat(Gc,S) 

Now, let us describe another approach to Chern-Simons gauge theory, based on quantiza- 
tion of moduli space spaces of flat connections. We focus on a theory with complex gauge 
group, which is the main subject of the present paper. ^ 

As we already mentioned in the introduction, Chern-Simons gauge theory (with any 
gauge group) associates a Hilbert space TCy: to a closed Riemann surface S and a vector in 
Tij] to every 3-manifold M with boundary S. We denote this vector as \M) G ?{■£. If there 
are two such manifolds, M-(_ and M_, glued along a common boundary S (with matching 
orientation), then the quantum invariant Z{M) that Chern-Simons theory associates to 
the closed 3-manifold M = Mj^ Us Af_ is given by the inner product of two vectors \Mj^) 
and |M_) in Hy, 

Z{M) = {M+\M^). (2.19) 

Therefore, in what follows, our goal will be to understand Chern-Simons gauge theory on 
manifolds with boundary, from which invariants of closed manifolds without boundary can 
be obtained via ( |2.19| ). 

Since the Chern-Simons action (1.1) is first order in derivatives, the Hilbert space "He 



is obtained by quantizing the classical phase space, which is the space of classical solutions 
on the 3-manifold M x S. According to ( |1.3| ), such classical solutions are given precisely 
by the flat connections on the Riemann surface S. Therefore, in a theory with complex 
gauge group Gc, the classical phase space is the moduli space of flat Gc connections on S, 
modulo gauge equivalence, 

A^flat(Gc,S) = Hom(7ri(S),Gc)/conj. (2.20) 



We will see, however, that quantization of Chern-Simons theory with complex gauge group Gc is closely 



related to quantization for compact group G. Indeed, as we shall describe in Section 



2.2.1 



this relation 



may be used to justify the "analytic continuation" approach to computing perturbative Gc invariants. 
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As a classical phase space, A^flatCGcS) comes equipped with a symplectic structure uo, 



which can also be deduced from the classical Chern-Simons action (1.1). Since we are inter- 
ested only in the "holomorphic" sector of the theory, we shall look only at the kinetic term 
for the field A (and not A); it leads to a holomorphic symplectic 2-form on A4flat(G'c, S): 

CO = / Tt SAaSA. (2.21) 

We note that this symplectic structure does not depend on the complex structure of S, 
in accord with the topological nature of the theory. Then, in Chern-Simons theory with 
complex gauge group Gc, the Hilbert space TCy: is obtained by quantizing the moduli space 
of flat Gc connections on S with symplectic structure ( |2.21 ): 



quantization of (A4flat(Gc,S),u;) -w Ht. ■ (2.22) 

Now, let us consider a closed 3-manifold M with boundary T, = dM, and its associated 
state \M) G "Hs- In a (semi-)classical theory, quantum states correspond to Lagrangian 
submanifolds of the classical phase space. Recall that, by definition, a Lagrangian subman- 
ifold £ is a middle-dimensional submanifold such that the restriction of to £ vanishes, 

u;\c = ■ (2.23) 

For the problem at hand, the phase space is 7Wflat(G'c, S) and the Lagrangian submanifold 
C associated to a 3-manifold M with boundary S = dM consists of the classical solutions 
on M. Since the space of classical solutions on M is the moduli space of flat Gc connections 
on M, 

A^flat(Gc, M) = Hom(7ri(M), Gc)/conj. , (2.24) 

it follows that 

£ = i{Ma.t{Gc,M)) (2.25) 
is the image of A^flat(G'C) -^) under the map 

i:A^fiat(Gc,M) ^A^flat(Gc,S) (2.26) 

induced by the natural inclusion map 7ri(S) — > 7ri(M). One can show that £ C A^flat(Gc, S) 
is indeed Lagrangian with respect to the symplectic structure ( |2.21| ). 

Much of what we described so far is very general and has an obvious analogue in Chern- 
Simons theory with arbitrary gauge group. However, quantization of Chern-Simons theory 
with complex gauge group has a number of good properties. In this case the classical phase 
space A^fiat(G£,S) is an algebraic variety; it admits a complete hyper-Kahler metric p8[ , 
and the Lagrangian submanifold £ is a holomorphic subvariety of A^flat(Gc, S). The 
hyper-Kahler structure on A^flat(Gc,S) can be obtained by interpreting it as the moduli 
space M.h{G, S) of solutions to Hitchin's equations on S. Note that this requires a choice 
of complex structure on S, whereas A1fiat(Gc,5]) = M.h{G,T,) as a complex symplectic 
manifold does not. Existence of a hyper-Kahler structure on A^flat(Gc,5]) considerably 
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simplifies the quantization problem in any of the existing frameworks, such as geometric 



quantization |54], deformation quantization or "brane quantization" |24]. 

The hyper-Kahler moduli space S) has three complex structures that we de- 

note as /, J, and K = IJ, and three corresponding Kahler forms, loj, ujj, and ujk- In the 
complex structure usually denoted by J, M.h{G, S) can be identified with TWflatCGc, 5^) as 
a complex symplectic manifold with the holomorphic symplectic form ( ^.21| ), 

oj = \{ujK + wi)- (2.27) 

Moreover, in this complex structure, C is an algebraic subvariety of A^flat(G'c, S). To be 
more precise, it is a (finite) union of algebraic subvarieties, each of which is defined by 
polynomial equations Ai = 0. In the quantum theory, these equations are replaced by 
corresponding operators Ai acting on Ti^, that annihilate the state \M). 

Now, let us consider in more detail the simple but important case when S is of genus 
1, that is S = T^. In this case, '7ri(S) = Z x Z is abelian, and 

A^flat(G'c,T2) = (TcxTc)/W, (2.28) 

where Tc is the maximal torus of Gc and W is the Weyl group. We parametrize each copy 
of Tc by complex variables / = (/i, . . . , Z^) and m = (mi, . . . , m^), respectively. Here, r is 
the rank of the gauge group Gc- The values of / and m are eigenvalues of the holonomies 
of the flat Gc connection over the two basic 1-cycles of S = T^. They are defined up to 
Weyl transformations, which act diagonally on Tc x Tc- 

The moduli space of fiat Gc connections on a 3-manifold M with a single toral bound- 
ary, dM = T^, defines a complex Lagrangian submanifold 

£ C (Tc X Tc)/W . (2.29) 

(More precisely, this Lagrangian submanifold comprises the top-dimensional (stable) com- 
ponents of the moduli space of flat Gc connections on S.) In particular, a generic irreducible 
component of C is defined by r polynomial equations 

Ai{l,m) = 0, i = 1,... ,r, (2.30) 

which must be invariant under the action of the Weyl group W (which simultaneously acts 
on the eigenvalues li, . . . ,1^ and mi, . . . , m^). In the quantum theory, these equations are 



replaced by the operator equations ( 1.12 ) 



Ai{l,rh) Z{M) = 0, i = l,...,r. (2.31) 
For S = the complex symplectic structure ( |2.21 ) takes a very simple form 



r 

I 



uj = -^dui Advi , (2.32) 



h 

i=l 
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where we introduce new variables u and v (defined modulo elements of the cocharacter 
lattice Acochar = Hom(C/(l), T)), such that / = —e" and m = e". In the quantum theory, 
u and V are replaced by operators u and v that obey the canonical commutation relations 

[ui,Vj\ = -h6ij . (2.33) 

As one usually does in quantum mechanics, we can introduce a complete set of states \u) 
on which u acts by multiplication, Ui\u) = Ui\u). Similarly, we let \v) be a complete basis, 
such that Vi\v) = Vi\v). Then, we can define the wave function associated to a 3-manifold 
M either in the n-space or u-space representation, respectively, as {u\M) or {v\M). We 
will mostly work with the former and let Z(M; u) := {u\M). 

We note that a generic value of u does not uniquely specify a flat Gc connection on M 
or, equivalently, a unique point on the representation variety ( 2.25| ). Indeed, for a generic 



value of u, equations ( 2.3C| ) may have several solutions that we label by a discrete parameter 



a. Therefore, in the u-space representation, flat Gc connections on M (previously labeled 
by the homomorphism p £ C) are now labeled by a set of continuous parameters u = 
(ui, . . . ,Ur) and a discrete parameter a: 

p i — > {u,a) . (2.34) 



The perturbative Gc invariant Z^p^ (M; h) can then be written in this notation as Z^'^^ (M; h, u) 
{u,a\M). Similarly, the coefficients in the /i-expansion can be written as sii'\u). 

To summarize, in the approach based on quantization of 7Wflat(Gc,S) the calcula- 
tion of Z^P\M;h) reduces to two main steps: i) the construction of quantum operators 
Ai{i,m), and ii) the solution of Schrodinger-like equations ( |2.31| ). Below we explain how 
to implement each of these steps. 

2.2.1 Analytic continuation 

For a generic 3-manifold M with boundary S = T^, constructing the quantum operators 
Ai{l, m) may be a difficult task. However, when M is the complement of a knot K in the 
3-sphere, 

M = S^\K, (2.35) 



there is a simple way to find the vlj's. Indeed, according to ( [1.17 ) these operators annihilate 



the polynomial knot invariants Pg,r{K; q), which are defined in terms of Chern-Simons 
theory with compact gauge group G, 

Aiii,m) PG,R{K;q)=0. (2.36) 

The operator /j acts on the set of polynomial invariants {Pg,r{K; q)} by shifting the highest 
weight A = (Ai, . . . , A^) of the representation R by the i-th basis elements of the weight 
lattice Awt, while the operator rhj acts simply as multiplication by q^^^'^. Let us briefly 
explain how this comes about. 

In general, the moduli space A^fiat(Gc,S) is a complexification of A^flat(G, S). The 
latter is the classical phase space in Chern-Simons theory with compact gauge group G and 



-17- 



can be obtained from A^flatCGc S) by requiring all the holonomies to be "real," i.e. in G. 
Similarly, restricting to real holonomies in the definition of C produces a Lagrangian sub- 
manifold in Ai{iat{G, S) that corresponds to a quantum state \M) in Chern-Simons theory 
with compact gauge group G. In the present example of knot complements, restricting to 
such "real" holonomies means replacing Tc by T in ( |2.28| ) and taking purely imaginary 
values of Ui and Vi in equations ( p.30D . Apart from this, the quantization problem is es- 
sentially the same for gauge groups G and Gc- In particular, the symplectic structure 
( |2.32 ) has the same form (with imaginary Ui and Vi in the theory with gauge group G) and 



the quantum operators Ai{l,'m) annihilate both PG,R{K;q) and Z^°'\M;h,u) computed, 
respectively, by Chern-Simons theories with gauge groups G and Gc- 

In order to understand the precise relation between the parameters in these theories, let 
us consider a Wilson loop operator, Wii{K), supported on iiT C in Chern-Simons theory 
with compact gauge group G. It is labeled by a representation R = Rx of the gauge group 
G, which we assume to be an irreducible representation with highest weight A € t^. As we 
already mentioned earlier (c/. eq. ( |1.13| )), the path integral in Chern-Simons theory on 
with a Wilson loop operator Wr{K) computes the polynomial knot invariant PcniK; q), 



with q = e . Using (2.19|), we can represent this path integral as {R\M), where \R) is the 



result of the path integral on a solid torus containing a Wilson loop Wji{K), and \M) is 
the path integral on its complement, M. 

In the semi-classical limit, the state \R) corresponds to a Lagrangian submanifold of 
A1flat(G, T^) = (T X T)/W defined by the fixed value of the holonomy m = on a small 
loop around the knot, where u is given by ( |1.14| ). The relation between m = e", which is an 
element of the maximal torus T of G, and the representation Rx is given by the invariant 
quadratic form — Tr (restricted to t). Specifically, 

m = exp(/iA*) G T, (2.37) 

where A* is the unique element of t such that A*(x) = — TrAx for all x £ t. In ( |1.15| ), we 
analytically continue this relation to m G Tc. 

For a given value of m = e", equations ( |2.3[1| ) have a finite set of solutions /„, labeled 
by a. Only for one particular value of a is the perturbative Gc invariant Z^'^\M;h,u) 
related to the asymptotic behavior of Pg,r{K; q). This is the value of a which maximizes 
Re {S^o"'\M;u)). For this a, we have (p^ ): 

— „ — = asymptotic expansion of Pg R(iv; . (2.38) 

For hyperbolic knots and u sufficiently close to 0, this "maximal" value of a is always 
a = geom. 

It should nevertheless be noted that the analytic continuation described here is not 
as "analytic" as it sounds. In particular, the limit ( 1.14| ) is very subtle and requires much 



care. As explained in [^], in taking this limit it is important that values of g = e"^^ avoid 
roots of unity. If one takes the limit with G ^rZ, which corresponds to the allowed 
values of the coupling constant in Chern-Simons theory with compact gauge group G, then 
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one can never see the exponential asymptotics ( [L.8| ) with Im(S'Q''^) > 0. The exponential 
growth characteristic to Chern-Simons theory with complex gauge group emerges only in 
the limit with h = u/\* and u generic. 

2.2.2 A hierarchy of differential equations 

The system of Schrodinger-like equations ( 2.31| ) determines the perturbative Gc invariant 
u) up to multiplication by an overall function of h, which can be fixed by suitable 
boundary conditions. 

In order to see in detail how the perturbative coefficients sii'\u) may be calculated 
and to avoid cluttering, let us assume that r = 1. (A generalization to arbitrary values of 
r is straightforward.) In this case, A{l,m) is the so-called A-polynomial of M, originally 
introduced in , and the system (|2.31| ) consists of a single equation 



A{1, m) z'-"'\M;h, u) 



0. 



(2.39) 



In the u-space representation the operator m = exp(u) acts on functions of u simply via 

d 



multiplication by e", whereas I = exp(?) + ivr) = exp(/i ^) acts as a "shift operator": 



^f{u) = e«/(n), lf{u) = f{u + h). 
In particular, the operators I and m obey the relation 

Im = q2ml, 

which follows directly from the commutation relation ( 2.33| ) for u and v, with 

„ _ „2h 



(2.40) 



(2.41) 



(2.42) 



We would like to recast eq. ( 2.39 ) as an infinite hierarchy of differential equations that 
can be solved recursively for the perturbative coefficients sii'\u). Just like its classical 
limit A{l,m), the operator A{l,7h) is a polynomial in /. Therefore, pushing all operators / 
to the right, we can write it as 



A{l,m) = 2^aj{m,h)P (2.43) 
i=o 

for some functions aj{m,h) and some integer d. Using ( 2.4C| ), we can write eq. ( 2.39| ) as 

(2.44) 

Then, substituting the general form (2.1) of Z^^^M; h,u), we obtain the equation 



^aj{m,h)Z^'^\M;h,u + jh) = 0. 



aj{m, h) exp 

i=o 



1 1 °° 



n=0 



. 



(2.45) 
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Since <5^°^ is independent of u, we can just factor out the —^(5^"^ log h term and remove it 
from the exponent. Now we expand everything in h. Let 

oo 

aj{m,h) = ^aj^p{m)h^ (2.46) 

p=0 



and 



oo 



-1 m= 



' _^ {r — m] 



suppressing the index a to simphfy notation. We can substitute p^ ) and ( p^ ) into 



(|2.45D and divide the entire expression by exp( K^Sn+iiu) ). The hierarchy of equations 
then follows by expanding the exponential in the resulting expression as a series in h and 
requiring that the coefficient of every term in this series vanishes. The first four equations 
are shown in Table |l[ 



iJ^'^ + ^^r + y^^'+J^U ) = 



Table 1: Hierarchy of differential equations derived from m) Z'^"^(M; ft, u) — 0. 

The equations in Table || can be solved recursively for the S'„(n)'s, since each Sn first 
appears in the (n + 1)'^* equation, differentiated only once. Indeed, after So is obtained, the 
remaining equations feature the S„>i linearly the first time they occur, and so determine 
these coefficients uniquely up to an additive constant of integration. 

The first equation, however, is somewhat special. Since ajfi{m) is precisely the coeffi- 
cient of P in the classical A-polynomial A{l,m), we can rewrite this equation as 

^(e^o(«),e") = 0. (2.48) 

This is exactly the classical constraint A{e'"^'^'^ , e") = that defines the complex Lagrangian 
submanifold £, with S'q{u) = v + iir. Therefore, WG CQiii integrate ctlong a, branch. (^Iq, — 
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e'"°'~^^'^,m = e") of £ to get the value of the classical Chern-Simons action ( |1.9D , 

= const + y 9\c, (2.49) 

where 9\c denotes a restriction to C of the Liouville 1-form on 7Wflat(G'c, S), 

e = vdu + indu . (2.50) 



The expression ( 2.49 ) is precisely the semi-classical approximation to the wave function 
Z^"\M;h,u) supported on the Lagrangian submanifold C, obtained in the WKB quan- 
tization of the classical phase space A^flat(G'c, S). By definition, the Liouville form 6 
(associated with a symplectic structure uj) obeys dO = ifiuj, and it is easy to check that this 
is indeed the case for the forms uj and 9 on Alflat(Gc,S) given by eqs. ( 2.32| ) and ( 2.5C| ), 
respectively. 



The semi-classical expression ( 2.49| ) gives the value of the classical Chern-Simons func- 



tional (^]^) evaluated on a flat gauge connection A^'^^ labeled by a homomorphism p. As 
we explained in ( |2.34| ), the dependence on p is encoded in the dependence on a continu- 
ous holonomy parameter u, as well as a discrete parameter a that labels different solutions 
Voi{u) to A[e'"'^™ , e") = 0, at a fixed value of u. In other words, a labels different branches of 
the Riemann surface A{1, m) = 0, regarded as a cover of the complex plane C parametrized 
by m = e", 

A(/o,m)=0. (2.51) 

Since A{l,m) is a polynomial in both I and m, the set of values of a is finite (in fact, its 
cardinality is equal to the degree of A{l,m) in I). Note, however, that for a given choice of 
a there are infinitely many ways to lift a solution la{u) to Va{u)] namely, one can add to 
Va{u) any integer multiple of 27ri. This ambiguity implies that the integral (|2.49| ) is defined 
only up to integer multiples of 27rm, 

/u pu 
log la{u ) du = / Va{u ) du + iTTU (mod 27rm) . (2.52) 

In practice, this ambiguity can always be fixed by imposing suitable boundary conditions 
on Sq"\u), and it never affects the higher-order terms sii'\u). Therefore, since our main 
goal is to solve the quantum theory (to all orders in perturbation theory) we shall not worry 



about this ambiguity in the classical term. As we illustrate later (see Section 4.2.1 ), it will 
always be easy to fix this ambiguity in concrete examples. 

Before we proceed, let us remark that if M is a hyperbolic 3-manifold with a single 
torus boundary E = dM and ^(s*^"™) is the "geometric" flat SL(2, C) connection associated 
with a hyperbolic metric on M (not necessarily geodesically complete), then the integral 
(2.52) is essentially the complexified volume function, i(Vol (M;u) + iCS(M;u)), which 



combines the (real) hyperbolic volume and Chern-Simons invariants of M. Specifically, 



*For example, imaginary it parametrizes a conical singularity. See, e.g. [ p6| , ^ and our discussion 
in Section B.l for more detailed descriptions of Vol (Af ; li) and CS(Af;u). In part of the literature {e.g. in 
pil), the parameters (u, v) are related to those used here by 2Mhere = Mthere and 2(i;hero + itt) = iithcro- We 
include a shift of ivr in our definition of v so that the complete hyperbolic structure arises at {u, v) — (0, 0). 
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the relation is |2|, |2|: 



= \ Vol (M; u) + iCS(M; u) + v^^oraiu) Re (u) + mu , (2.53) 
modulo the integration constant and multiples of Iniu. 
2.3 Classical and quantum symmetries 

A large supply of 3-manifolds with a single toral boundary can be obtained by considering 
knot complements ( |2.35| ); our main examples in Section Q are of this type. As discussed 
above, the Lagrangian subvariety C £ A^flat(G'c, S) for any knot complement M is defined 
by polynomial equations ( 2.3C| ). Such an C contains multiple branches, indexed by a, 



corresponding to the different solutions to {Ai{l,m) = 0} for fixed m. In this section, we 
describe relationships among these branches and the corresponding perturbative invariants 
using the symmetries of Chern-Simons theory with complex gauge group 

Before we begin, it is useful to summarize what we already know about the branches 
of £. As mentioned in the previous discussion, there always exists a geometric branch — 
or in the case of rank r > 1 several geometric branches — when M is a hyperbolic knot 
complement. Like the geometric branch, most other branches of C correspond to genuinely 
nonabelian representations p : tti{T,) Gc. However, for any knot complement M there 
also exists an "abelian" component of C, described by the equations 

1-^ = ... = l^ = l. (2.54) 



Indeed, since Hi{M) is the abelianization of 7ri(M), the representation variety ( p. 24 ) always 



has a component corresponding to abelian representations that factor through Hi{M) = Z, 

7ri(M) ^ Hi{M) Gc ■ (2.55) 

The corresponding flat connection, ^(^'^^^) , is characterized by the trivial holonomy around 
a 1-cycle of S = which is trivial in homology Hi (M) ; choosing it to be the 1-cycle whose 
holonomy was denoted by / = {h, . . . ,lr) we obtain ( 2.54| ). Note that, under projection to 



the u-space, the abelian component of £. corresponds to a single branch that we denote by 
a = abel. 

The first relevant symmetry of Chern-Simons theory with complex gauge group Gc is 
conjugation. We observe that for every flat connection A^^^ on M, with p = (u, a), there 
is a conjugate flat connection A^^^ := A^P^ corresponding to a homomorphism p = {u,a). 
We use a to denote the branch of C "conjugate" to branch a; the fact that branches 



of C come in conjugate pairs is reflected in the fact that eqs. ( 2.30 ) have real (in fact, 
integer) coefficients. The perturbative expansions around A^ and A^ are very simply 
related. Namely, by directly conjugating the perturbative path integral and noting that the 
Chern-Simons action has real coefficients, we find^ Z'^{M;h,u) = Z^"\M;h,u). The latter 
^ More explicitly, letting Icsih,A) = -jj-^Jj^jTr[AAA+lAAAAA),we have 

Z-(M;ft,u) = ( / VAe''^''^'''-^A =f VAe'^^^'''-^^ = Z^''\M;h,u) . 

I (u,a) J J (u,q) 
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partition function is actually in the antiholomorphic sector of the Chern-Simons theory, but 
we can just rename {h, u) i— > (h, u) (and use analyticity) to obtain a perturbative partition 
function for the conjugate branch in the holomorphic sector, 



Z^^\M;n,u) = Z(.'^){M;h,u) . (2.56) 



Here, for any function f{z) we define f{z) := f{z). In particular, if / is analytic, f{z) = 
^ /„2:", then / denotes a similar function with conjugate coefficients, f{z) = ^ fnZ^- 

In the case r = 1, the symmetry ( ^.56| ) implies that branches of the classical SL{2,C) 
A-polynomial come in conjugate pairs Va and Va{u) = v^{u). Again, these pairs arise 
algebraically because the A-polynomial has integer coefficients. (See e.g. |^ § for a 
detailed discussion of properties of A{l,m).) Some branches, like the abelian branch, may 
be self-conjugate. For the abelian branch, this is consistent with ^^^^^'^ = q. The geometric 



branch, on the other hand, has a distinct conjugate because Vol(M;0) > 0; from ( 2.53| ) 
we see that its leading perturbative coefficient obeys 

5f °''j)(M) = ^ - Vol (M; u) + iCS(M; u) + v^{u)Re (u) - iiru . (2.57) 
In general, we have 



5f)(n) = 5J"^(m) (mod27rn). (2.58) 

Now, let us consider symmetries that originate from geometry, i.e. symmetries that 
involve involutions of M, 

r: M^M. (2.59) 
Every such involution restricts to a self-map of S = dM, 

t|s: S^S, (2.60) 

which, in turn, induces an endomorphism on homology, HiiTi). Specifically, let us con- 
sider an orientation-preserving involution r which induces an endomorphism (—1,-1) on 
Hi{T?) = Z X Z. This involution is a homeomorphism of M; it changes our definition of 
the holonomies, 

rrii — and li ^ — , (2.61) 

rrii k 

leaving the symplectic form ( p.32| ) invariant. Therefore, it preserves both the symplectic 
phase space A^fiat(G'c, S) and the Lagrangian submanifold C (possibly permuting some of 
its branches). 

In the basic case of rank r = 1, the symmetry ( p. 61 ) corresponds to the simple, well- 



known relation A(l~^ ,m~^) = A{l,m), up to overall powers of I and m. Similarly, at 
the quantum level, A{1~^ ,rh^^) = A{i,m) when A{i,m) is properly normalized. Branches 
of the A-polynomial are individually preserved, implying that the perturbative partition 
functions (and the coefficients S^^) are all even: 

z'^°'\M;h,-u) = z'^''\M;h,u), (2.62) 
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modulo factors of e^'^^l^ that are related to the ambiguity in Sq°'\u). Note that in the 
r = 1 case one can also think of the symmetry (2.61) as the Weyl reflection. Since, by 



definition, holonomies that differ by an element of the Weyl group define the same point in 
the moduli space ( 2.28^ , it is clear that both A4fiat(G'c, S) and £ are manifestly invariant 



under this symmetry. (For r > 1, Weyl transformations on the variables / and m lead to 
new, independent relations among the branches of £.) 

Finally, let us consider a more interesting "parity" symmetry, an orientation-reversing 
involution 

P : M (2.63) 

P|e : S^S (2.64) 

that induces a map (1,-1) on -ffi(S) = Z x Z. This operation by itself cannot be a 
symmetry of the theory because it does not preserve the symplectic form ( |2.32| ). We can 
try, however, to combine it with the transformation h — > —h to get a symmetry of the 
symplectic phase space (A^flat(Gc, S), a;). We are still not done because this combined 
operation changes the orientation of both S and M, and unless M = M the state |M) 
assigned to M will be mapped to a different state \M). But if M is an amphicheiral^ 
manifold, then both "T-i — *■ ^ and k ^ j- (independently) become symmetries of the 
theory, once combined with h — > —h. This now implies that solutions come in signed pairs, 
Va{u) and fa(n) = —Vaiu), such that the corresponding perturbative Gc invariants satisfy 

(M; h, u) = (M; -h,u) . (2.65) 

For the perturbative coefficients, this leads to the relations 

S^^\u) = (mod 2Tru) , (2.66) 

Si''\u) = i-ir+'si-\u) n>l. (2.67) 

Assuming that the ~ i|Vol (M; u) behavior of the geometric and conjugate branches 
is unique, their signed and conjugate pairs must coincide for amphicheiral 3-manifolds. 
Sll''°'^\ then, is an even analytic function of u with strictly real series coefficients; at 
u £ iM, the Chern-Simons invariant CS{M;u) will vanish. 

2.4 Brane quantization 

Now, let us briefly describe how the problem of quantizing the moduli space of flat connec- 
tions, A^flat(Gc, S), would look in the new approach [24| based on the topological A-model 



and D-branes. Although it can be useful for a better understanding of Chern-Simons the- 
ory with complex gauge group, this discussion is not crucial for the rest of the paper and 
the reader not interested in this approach may skip directly to Section 0. 



In the approach of |24], the problem of quantizing a symplectic manifold N with 
symplectic structure u is solved by complexifying (iV, uj) into {Y, Q) and studying the A- 
model of Y with symplectic structure = ImQ. Here, y is a complexification of N, i.e. 



®A manifold is called chiral or amphicheiral according to whether the orientation cannot or can be 
reversed by a self-map. 
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a complex manifold with complex structure X and an antiholomorphic involution 



t:Y^Y, (2.68) 

such that is contained in the fixed point set of r and t*T = — T. The 2-form on y is 
holomorphic in complex structure X and obeys 

T*J] = n (2.69) 

and 

n|7v = a;. (2.70) 

In addition, one needs to pick a unitary line bundle £ ^ y (extending the "prequantum line 
bundle" £ — > A^) with a connection of curvature Reil. This choice needs to be consistent 
with the action of the involution r, meaning that r : y — > y lifts to an action on £, such 
that tItv = id. To summarize, in brane quantization the starting point involves the choice 
of y, fi, £, and r. 

In our problem, the space N = A4fiat(Gc,S) that we wish to quantize is already 
a complex manifold. Indeed, as we noted earlier, it comes equipped with the complex 
structure J (that does not depend on the complex structure on E). Therefore, its com- 
plexification^ is y = x with the complex structure on A'^ being prescribed by — J and 
the complex structure on Y being X = (J, — J). The tangent bundle TY = TN © TN is 
identified with the complexified tangent bundle of A^, which has the usual decomposition 
CTN = T^^^N © T^^^N. Then, the "real shoe" A^ is embedded in Y as the diagonal 

N ^ {x,x) e N yil^ . (2.71) 

In particular, N is the fixed point set of the antiholomorphic involution t -.Y ^Y which 
acts on {x,y) &Y as t : {x^y) ^-^ {y,x). 

Our next goal is to describe the holomorphic 2-form O that obeys ( 2.69] ) and ( 2.7C| ) 
with^ 

uj = —{iVK + iuJi) - —{ujK - iioi) (2.72) 

Note, that {ujk + i^i) is holomorphic on A^ = 7Wfiat(Gc, S) and {uk — i^i) is holomorphic 
on A^. Moreover, if we take t to be a complex conjugate of t, the antiholomorphic involution 
r maps t{ujK + i^i) to i{ujK — i^i)-, so that t*uj = w. Therefore, we can simply take 

O = ^(a;J) + - ^^{J^ - ^u:f ) (2.73) 

where the superscript i = 1,2 refers to the first (resp. second) factor in y = A^ x A^. It 
is easy to verify that the 2-form J7 defined in this way indeed obeys (Imr2)~^Rer2 = X. 



^Notice, since in our problem we start with a hyper-Kahler manifold A^flat(Gc, S), its complexification Y 
admits many complex structures. In fact, Y has holonomy group Sp{n) x Sp{n), where n is the quaternionic 
dimension of A^flat(G'c, S). 

^Notice, while in the rest of the paper we consider only the "holomorphic" sector of the theory (which is 
sufficient in the perturbative approach), here we write the complete symplectic form on A1flat(Gc,S) that 



follows from the classical Chern-Simons action (1.1), including the contributions of both fields A and A. 
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Moreover, one can also check that if i is a complex conjugate of t then the restriction of 
COY = Imfi to the diagonal ( |2.71| ) vanishes, so that the "real slice" C y, as expected, is 
a Lagrangian submanifold in (Y,ujy)- 

Now, the quantization problem can be realized in the A-model of Y with symplectic 
structure ujy = ImO. In particular, the Hilbert space Tis is obtained as the space of 
{Bcc,B') strings, 

We = space of {Bcc,B') strings (2.74) 

where Bcc and B' are A-branes on Y (with respect to the symplectic structure wy). The 
brane B' is the ordinary Lagrangian brane supported on the "real slice" N CY. The other 
A-brane, Bcc, is the so-called canonical coisotropic brane supported on all of Y. It carries 
a Chan-Paton line bundle of curvature F = ReO. Note that for [F] to be an integral 
cohomology class we need Re (t) £ Z. Since in the present case the involution r fixes the 
"real slice" pointwise, it defines a hermitian inner product on TC^, which is positive definite. 



3. A state integral model for perturbative SL{2, C) Chern-Simons theory 

In this section, we introduce a "state integral" model for Z^^^ (M; h) in the simplest case 



of Gc = SL{2,C). Our construction will rely heavily on the work of Hikami pq , 27 1, 
where he introduced an invariant of hyperbolic 3-manifolds using ideal triangulations. The 
resulting invariant is very close to the state integral model we are looking for. However, in 
order to make it into a useful tool for computing Z^p\M;fi) we will need to understand 
Hikami's construction better and make a number of important modifications. In particular, 
as we explain below, Hikami's invariant is defined as a certain integral along a path in the 
complex plane (or, more generally, over a hypersurface in complex space) which was not 
specified^ in the original work (2^, Another issue that we need to address is how to 
incorporate in Hikami's construction a choice of the homomorphisms ( |l.4D , 



p : 7ri(M) ^ 5L(2,C). (3.1) 

(The original construction assumes very special choices of p that we called "geometric" 
in Section §.) It turns out that these two questions are not unrelated and can be ad- 
dressed simultaneously, so that Hikami's invariant can be extended to a state sum model 
for Z^p\M] h) with an arbitrary p. 

Throughout this section, we work in the u-space representation. In particular, we use 
the identification ( |2.34| ) and denote the perturbative SL(2,C) invariant as Z^°'\M]h,u). 

3.1 Ideal triangulations of hyperbolic 3-manifolds 

The construction of a state integral model described in the rest of this section applies to 
orientable hyperbolic 3-manifolds of finite volume (possibly with boundary) and uses ideal 



®One choice, briefly mentioned in p6| |27| , is to integrate over the real axis (resp. real subspace) of the 
complex parameter space. While this choice is in some sense natural, a closer look shows that it cannot be 
the right one. 
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triangulations in a crucial way. Therefore, we begin this section by reviewing some relevant 
facts from hyperbolic geometry (more details can be found in | 



47, 4 



Recall that hyperbolic 3-space can be represented as the upper half-space {{xi, X2, xs) 
X3 > 0} with metric ( |2.7| ) of constant curvature —1. The boundary dM.^, topologically an 
S^, consists of the plane X3 = together with the point at infinity. The group of isometrics 
of H'^ is PSL{2,C), which acts on the boundary via the usual Mobius transformations. In 
this picture, geodesic surfaces are spheres of any radius which intersect dM^ orthogonally. 

An ideal tetrahedron A in has by definition all its faces along geodesic surfaces, and 
all its vertices in dM.^ — such vertices are called ideal points. After Mobius transformations, 
one can fix three of the vertices at (0,0,0), (1,0,0), and infinity. The coordinate of the 
fourth vertex {xi,X2,0), with X2 > 0, defines a complex number z = xi + ix2 called the 
shape parameter (sometimes also called edge parameter). At various edges, the faces of the 
tetrahedron A form dihedral angles argZj (j = 1,2,3) as indicated in Figure |2|, with 

(3.2) 



Z2 



1 



Z3 



1 



The ideal tetrahedron is noncompact, but has finite volume given by 

Vol (A,) = D{z), 



(3.3) 



where D{z) is the Bloch-Wigner dilogarithm function, related to the usual dilogarithm Li2 
(see Section |3.3| ) by 

D{z) = Im (Li2(z)) + arg(l - z) log |z| . (3.4) 

Note that any of the Zj can be taken to be the shape parameter of A, and that D{zj) = 
Vol (Ax) for each j. We will allow shape parameters to be any complex numbers in C — 
{0, 1}, noting that for 2; G M an ideal tetrahedron is degenerate and that for Im z < it 
technically has negative volume due to its orientation. 



(projection to boundary plane) 





Figure 2: An ideal tetrahedron in H"^. 

A hyperbolic structure on a 3-manifold is a metric that is locally isometric to H^. A 3- 
manifold is called hyperbolic if it admits a hyperbolic structure that is geodesically complete 
and has finite volume. Most 3-manifolds are hyperbolic, including the vast majority of knot 
and link complements in S^. Specifically, a knot complement is hyperbolic as long as the 
knot is not a torus or satellite knot [E^]. Every closed 3-manifold can be obtained via Dehn 
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surgery on a knot in S^, and for hyperbolic knots all but finitely many such surgeries yield 
hyperbolic manifolds |p6[| . 

By the Mostow rigidity theorem |]37| , the complete hyperbolic structure on a hy- 
perbolic manifold is unique. Therefore, geometric invariants like the hyperbolic volume 
are actually topological invariants. For the large class of hyperbolic knot complements in 
S^, the unique complete hyperbolic structure has a parabolic holonomy with unit eigenval- 
ues around the knot. In SL{2,C) Chern-Simons theory, this structure corresponds to the 
"geometric" flat connection with u = 0. As discussed in Section hyperbolic 

manifolds with complete hyperbolic structures can also be described as quotients H'^/F. 

Given a hyperbolic knot complement, one can deform the hyperbolic metric in such a 
way that the holonomy u is not zero. Such deformed metrics are unique in a neighborhood 
of n = 0, but they are not geodesically complete. For a discrete set of values of u, one can 
add in the "missing" geodesic, and the deformed metrics coincide with the unique complete 
hyperbolic structures on closed 3-manifolds obtained via appropriate Dehn surgeries on the 
knot in S^. For other values of u, the knot complement can be completed by adding either 
a circle or a single point, but the resulting hyperbolic metric will be singular. For 
example, if n G iM one adds a circle and the resulting metric has a conical singularity. 
These descriptions can easily be extended to link complements (i.e. multiple cusps), using 
multiple parameters Uk, one for each link component. 

Any orientable hyperbolic manifold M is homeomorphic to the interior of a compact 
3-manifold M with boundary consisting of finitely many tori. (The manifold M itself can 
also be thought of as the union of M with neighborhoods of the cusps, each of the latter 
being homeomorphic to x [0,oo).) All hyperbolic manifolds therefore arise as knot or 
link complements in closed 3-manifolds. Moreover, every hyperbolic manifold has an ideal 
triangulation, i.e a finite decomposition into (possibly degenerate^'') ideal tetrahedra; see 



e.g. H 



To reconstruct a hyperbolic 3-manifold M from its ideal triangulation {Ajj^j^, faces 
of tetrahedra are glued together in pairs. One must remember, however, that vertices of 
tetrahedra are not part of M, and that the combined boundaries of their neighborhoods 
in M are not spheres, but tori. (Thus, some intuition from simplicial triangulations no 
longer holds.) There always exists a triangulation of M whose edges can all be oriented in 
such a way that the boundary of every face (shared between two tetrahedra) has two edges 
oriented in the same direction (clockwise or counterclockwise) and one opposite. Then the 
vertices of each tetrahedron can be canonically labeled 0, 1, 2, 3 according to the number 
of edges entering the vertex, so that the tetrahedron can be identified in a unique way with 
one of the two numbered tetrahedra shown in Figure ^ of the next subsection. This at 
the same time orients the tetrahedron. The orientation of a given tetrahedron Aj may not 
agree with that of M; one defines = 1 if the orientations agree and = — 1 otherwise. 
The edges of each tetrahedron can then be given shape parameters (z ^2*^ Zg*^ ) , running 
counterclockwise around each vertex (viewed from outside the tetrahedron) if = 1 and 
clockwise if e,- = —1. 



"it is conjectured and widely believed that nondegenerate tetrahedra alone are always sufficient. 
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For a given M with cusps or conical singularities specified by holonomy parameters Uk, 

(i) 

the shape parameters of the tetrahedra Aj in its triangulation are fixed by two sets of 

(i) 

conditions. First, the product of the shape parameters Zj at every edge in the triangulation 
must be equal to 1, in order for the hyperbolic structures of adjacent tetrahedra to match. 
More precisely, the sum of some chosen branches of log zj*'* (equal to the standard branch 
if one is near the complete structure) should equal 27rz, so that the total dihedral angle 
at each edge is 2tt. Second, one can compute holonomy eigenvalues around each torus 
boundary in M as a product of Zj s by mapping out the neighborhood of each vertex in 
the triangulation in a so-called developing map, and following a procedure illustrated in, 
e.g., 1^]. There is one distinct vertex "inside" each boundary torus. One then requires 
that the eigenvalues of the holonomy around the kth boundary are equal to 6^"*=. These 
two conditions will be referred to, respectively, as consistency and cusp relations. 



Every hyperbolic 3-manifold has a well-defined class in the Bloch group |41]. This is a 
subgroup^^ of the quotient of the free Z- module Z[C — {0, 1}] by the relations 



\y] + 



[^1 

X 



.1 



1-r 



i-y' 



0. 



(3.5) 



This five-term or pentagon relation accounts for the fact that a polyhedron with five ideal 
vertices can be decomposed into ideal tetrahedra in multiple ways. The five ideal tetrahedra 
in this polyhedron (each obtained by deleting an ideal vertex) can be given the five shape 
parameters x, y, y/x,... appearing above. The signs of the different terms correspond to 
orientations. Geometrically, an instance of the five-term relation can be visualized as the 
2-3 Pachner move, illustrated in Figure 0. 





Figure 3: The 2-3 Pachner move. 



The class [M] of a hyperbolic 3-manifold M in the Bloch group can be computed 
by summing (with orientation) the shape parameters [z] of any ideal triangulation, but 
it is independent of the triangulation. Thus, hyperbolic invariants of 3-manifolds may be 
obtained by functions on the Bloch group — i.e. functions compatible with (|3.5D . For 
example, the Bloch- Wigner function (|3.4| ) satisfies 



D{x) - D{y) + D{!^) - D 



,1 



l-y- 



+ D( 



.1 



0, 



(3.6) 



Namely, the kernel of the map [2] 22: A (1 — z) G 



Az ' 



acting on this quotient module. 
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and the hyperbolic volume of a manifold M triangulated by ideal tetrahedra {Aj}^^ can 
be calculated as 

N 

Yo\{M)=^eiD{z^^). (3.7) 

i=l 

The symbols here could be removed if shape parameters were assigned to tetrahedra in a 
manner independent of orientation, noting that reversing the orientation of a tetrahedron 
corresponds to sending z ^ 1/ z and that D{1/ z) = —D{z). This is sometimes seen in the 
literature. 

The complexified volume i(Vol (M) + iCS(M)) is trickier to evaluate. For a hyperbolic 
manifold with a spin structure, corresponding to full SL{2,C) holonomies, this invariant 
is defined modulo 47r^. Here, we will outline a computation of the complexified volume 
modulo vr^, following |40|; for the complete invariant modulo 47r^, see j6^. To proceed, one 
must first make sure that the three shape parameters z(^=zP,zi'^ , and are specifically 
assigned to edges [f[*\'U2*^], [i^i*'*, ^3*^], and [^j*\w2*^] (respectively) in each tetrahedron Aj 
of an oriented triangulation of M, where [va\vj^^] denotes the edge going from numbered 
vertex Va^ to numbered vertex One also chooses logarithms {'w^\ w.^^ ) of the 

shape parameters such that 

A'=±z^\ e-^'=±z«, A'=±z^\ (3.8a) 
wf^ + w!^ ^ + = V i , (3.8b) 

and defines integers (gW,rW) by 

wf^ = Log(z(^)) + niq^^ , = -Log(l - z^'^) + niA'^ , (3.9) 

where Log denotes the principal branch of the logarithm, with a cut from to — cxd. For a 
consistent labeling of the triangulation, called a "combinatorial flattening," the sum of log- 
parameters w^j^ around every edge must vanish, and the (signed^'^) sum of log-parameters 
along the two paths generating TriiT"^) = 1? for any boundary (cusp) must equal twice 
the logarithm of the SL{2, C) holonomies around these paths. ^'^ The complexified volume 
is then given, modulo tt^, as 

N 

i(Vol (M) + iCS(M)) = ^eiL(z«;g»,r«)- J] {vkUk + it^uu) , (3.10) 

i=l cusps fc 

with 

1 -K^ar TT^ 

L(z; g, r) = Li2(z) + -(Log(z) + 7rig)(Log(l - z) + mr) + - — . (3.11) 



Signs arise from tetrahedron orientations and the sense in which a path winds around edges; see ^Oj , Def. 4.2. 
^^Explicitly, in the notation of Section 2.2 and above, the sum of log-parameters along the two paths in 



the neighborhood of the fcth cusp must equal 2uk and 2vk + 2-Ki, respectively. 
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The function L(z;q,r), a modified version of the Rogers dilogarithm, satisfies a five-term 
relation in an extended Bloch group that hfts ( |3.5| ) in a natural way to the space of log- 
parameters.^"^ 



3.2 Hikami's invariant 



We can now describe Hikami's geometric construction. Roughly speaking, to compute the 
invariant for a hyperbolic manifold M, one chooses an ideal triangulation of M, assigns an 
infinite-dimensional vector space V or V* to each tetrahedron face, and assigns a matrix 
element in V <SiV fSiV* 0V* to each tetrahedron. These matrix elements depend on a small 
parameter h, and in the classical — > limit they capture the hyperbolic structure of the 
tetrahedra. The invariant of M is obtained by taking inner products of matrix elements 
on every pair of identified faces (gluing the tetrahedra back together), subject to the cusp 
conditions described above in the classical limit. 

To describe the process in greater detail, we begin with an orientable hyperbolic man- 
ifold that has an oriented ideal triangulation {Aj}^^, and initially forget about the hyper- 
bolic structures of these tetrahedra. As discussed in Section 3T and indicated in Figure |^, 
each tetrahedron comes with one of two possible orientations of its edges, which induces 
an ordering of its vertices wj*^ (the subscript j here is not to be confused with the shape 
parameter subscript in (|3.2|)), an ordering of its faces, and orientations on each face. The 
latter can be indicated by inward or outward-pointing normal vectors. The faces (or their 
normal vectors) are labelled by p^j \ in correspondence with opposing vertices. The normal 
face- vectors of adjacent tetrahedra match up head-to-tail (and actually define an oriented 
dual decomposition) when tetrahedra are glued to form M. 





'1 P4 



Figure 4: Oriented tetrahedra, to which matrix elements (pi,p3|S|p2,P4) (left) 
(p2,P4|S"^|pi,P3) (right) are assigned. 



and 



Given such an oriented triangulation, one associates a vector space V to each inward- 
oriented face, and the dual space V* to each outward-oriented face. (Physicists should think 



^''The branch of Li2 in (3.11) is taken to be the standard one, with a cut from 1 to +oo. Note, however 

bi 

/- 



that we could also take care of ambiguities arising from the choice of dilogarithm branch by rewriting (3.11) 

t dt 



in terms of the function L{w) = 
because all the residues of t/{l 
function ^2(6™) + w log(l - e") 



This is a well-defined holomorphic function 



e ) are integer multiplies of 27ri, and it coincides with a branch of the 
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of these spaces as "Hilbert" spaces obtained by quantizing the theory on a manifold with 
boundary.) The elements of V are represented by complex- valued functions in one variable, 
with adjoints given by conjugation and inner products given by integration. Abusing the 
notation, but following the very natural set of conventions of we denote these 

complex variables by p^j \ which we used earlier to label the corresponding faces of the 
tetrahedra. As a result, to the boundary of every tetrahedron Aj one associates a vector 

(i) (i) (i) 

space V ^ V (Si V* V* , represented by functions of its four face labels, p\ , P2 , p\ , 
and ^ . 

To each tetrahedron one assigns a matrix element {pi,Pz\^\p2,Pa) or {p2,PA\^ ^^\pi,P'i) , 
depending on orientation as indicated in Figure |^. Here, the matrix S acts on functions 
fiPi,P2) €V^V as 

S = e«iP2/2ft$;i(pi+g2-P2), (3.12) 
where pif = pif and (jif = 2h-^f. The function <^fi is the quantum dilogarithm, to 
be described in the next subsection. Assuming that / = / dp\p){p\ = 1 and 
{p\q) = (the exact normalizations are not important for the final invariant), one 

obtains via Fourier transform 

(Pi,P3|S|p2,P4) = , , $R(p4-P3 + »7r + ft)e='''V 6 7, (3.13) 

(p.,P4|S-b„P3) = L_^,4(.UP3-..)--.^) (3,^4) 

^J-^mh ^h(P4 - P3 - «7r - ft) 

In the classical limit — > 0, the quantum dilogarithm has the asymptotic ^h{p) ~ 
^Li2(— e^). One therefore sees that the classical limits of the above matrix elements look 
very much like exponentials of ^ times the complexified hyperbolic volumes of tetrahe- 
dra. For example, the asymptotic of ( p.l3|) coincides with exp(L(z; •, ■)/{2h)) if we identify 
gP4-P3 with z and e~^^^ with 1/(1 — z). For building a quantum invariant, however, only 
half of the variables pj really "belong" to a single tetrahedron. Hikami's claim [26, 27 1 is 
that if we only identify a shape parameter 

zW=eP4-P3 , (3.15) 

for every tetrahedron Aj, the classical limit of the resulting quantum invariant will com- 
pletely reproduce the hyperbolic structure and complexified hyperbolic volume on M.^^ 

To finish calculating the invariant of M, one glues the tetrahedra back together and 
takes inner products in every pair V and V* corresponding to identified faces. This amounts 



to multiplying together all the matrix elements (|3.13| ) or ( p. 14 ) , identifying the p^^^ variables 



on identified faces (with matching head-to-tail normal vectors), and integrating over the 
2A^ remaining p's. To account for possible toral boundaries of M, however, one must revert 
back to the hyperbolic structure. This allows one to write the holonomy eigenvalues {e"'- } 
as products of shape parameters {z^^\ 1 — l/z^^\ 1/(1 ~ -2^*^)), and, using ( ^.15 ), to turn 



every cusp condition into a linear relation of the form '^p^s= 2uk- These relations are 



Eqn. (3.15) is a little different from the relation appearing in ji^, |2^], because our convention for 



assigning shape parameters to edges based on orientation differs from that of ^ 
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then inserted as delta functions in the inner product integral, enforcing global boundary 
conditions. In the end, noting that each matrix element ( 3.13D -( ^.14| ) also contains a delta 
function, one is left with N — boCE) nontrivial integrals, where 60 (5^) is the number of 
connected components of S = dM. For example, specializing to hyperbolic 3-manifolds 
with a single torus boundary S = T'^, the integration variables can be relabeled so that 
Hikami's invariant takes the form 

1 r ^ 

H{M- h,u) = ^^^^^ J llM9iiP,'^u) + EiiiTT + h)y' e/(P-2-,?i)/2/l _ , , dp^_^ 

(3.16) 

The gi are linear combinations of (pi, . . . ,pjv_i,2n) with integer coefficients, and / is a 
quadratic polynomial, also with integer coefficients for all terms involving p^'s or u. In the 
classical limit, this integral can naively be evaluated in a saddle point approximation, and 
Hikami's claim is that the saddle point relations coincide precisely with the consistency 
conditions for the triangulation on M. There is more to this story, however, as we will see 
in Section 3.4 . 

3.3 Quantum dilogarithm 

Since quantum dilogarithms play a key role here, we take a little time to discuss some of 
their most important properties. 

Somewhat confusingly, there are at least three distinct — though related — functions 
which have occurred in the literature under the name "quantum dilogarithm" : 

i) the function Li2(x;g) defined for x, q gC with \q\ < 1 by 

00 „ 

whose relation to the classical dilogarithm function Li2(a;) = — 5- is that 

n=i 

U2{x;e^'') ^L^2(x) ash^O; (3.18) 

a) the function (x; q)oo defined for < 1 and all x G C by 

00 

(x;g)oo = Hil-q'^x), (3.19) 

T-=0 

which is related to Li2(x;g) for |x| < 1 by 

ix;q)oo = exp(-Li2(x;g)) ; (3.20) 

and finally 

in) the function ^(z; r) defined for Re (r) > and 2|Re (z)] < 1 + Re (r) by 

N /^l /" e^^^ dx\ 
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(here M^^' denotes a path from — oo to oo along the real line but deformed to pass over 
the singularity at zero), which is related to {x;q)oo by 



$(z;r) 



e(z + T/2); e(r))^ 



(_e((.-l/2)/r);e(-l/T))^ 
(-e((z + l/2)/r);e(l/r))^ 
(_e(z-r/2);e(-T))^ 



if Im (t) > 0, 



if Im (r) < 0. 



(3.22) 



' oo 

(Here and in future we use the abbreviation e(x) = e^'^*^.) 
It is the third of these functions, in the normalization 



$fif^) 



; — ) 

27ri ITT 



(3.23) 



which occurs in our "state integral" and which we will take as our basic "quantum dilog- 
arithm," but all three functions play a role in the analysis, so we will describe the main 
properties of all three here. We give complete proofs, but only sketchily since none of this 
material is new. For further discussion and proofs, see, e.g., jl^, 14, 15, 51, 2C], and |57| 
(subsection II.l.D). 

1. The asymptotic formula ( ^.18 ) can be refined to the asymptotic expansion 



X h 



U,(.; ~ Li,(x) - - l„g(l -.)-_-+ OS^ + + 



(3.24) 



as /i — > with X fixed, in which the coefficient of hP~^ for n > 2 is the product of 
—2'^~^Bn/nl (here Bn is the nth Bernoulli number) with the negative-index polylogarithm 
Li2_n(2;) G Q[]^] • More generally, one has the asymptotic formula 



E2" ^Bn{X) , , 
— Li2_„(x) h 

n=0 



nl 



(3.25) 



as ^ ^ with A fixed, where Bn{t) denotes the nth Bernoulli polynomial. Both formulas 
are easy consequences of the Euler-Maclaurin summation formula. By combining ( 3.22| ), 
(3.20), and ( 3.25| ), one also obtains an asymptotic expansion 



^n{z + 2Xh) = exp ^ 



°° 2"-ii?„(l/2 + A) 1 

— n Li2-„(— e 



nl 



(3.26) 



(To derive this, note that in ( |3.22| ) {—e{z ± l/2)/r; e(ibl/r))oo ~ 1 to all orders in h as 



2. The function (x; q)oo and its reciprocal have the Taylor expansions 

oo 

T — 



{x;q)c 



°^ n(„-l) 

2^ ^7Z^ Q ' ^ 



n=0 



iQ)n 



{x;q) 



oo 



n=0 



iQ)n 



X 



(3.27) 



^^This is the unique polynomial satisfying Bn{t)dt = i", and is a monic polynomial of degree 

with constant term _B„. 
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around x = 0, where 



{Q)n 



(3.28) 



is the nth (j-Pochhammer symbol. These, as well as formula ( |3.2C1| ), can be proved easily 
from the recursion formula (x; q)oo = (1 — x){qx; q)oo, which together with the initial value 
{0',<i)oo = 1 determines the power series {x;q)oo uniquely. (Of course, ( ^.20 ) can also be 
proved directly by expanding each term in log(l — q'^x) as a power series in x.) Another 



famous result, easily deduced from ( 3.27| ) using the identity Yl 

m-n=k [QjmiQjn {Q)c 



for 



all /c G Z, is the Jacobi triple product formula 

oo 



(-l)'=g 2 x^ 



(3.29) 



relating the function (x; q)oo to the classical Jacobi theta function. 

3. The function ^{z;t) defined (initially for Re (r) > and |Re(z)| < ^ + |Re(r)) 



by ( 3.21 ) has several functional equations. Denote by I{z; r) the integral appearing in this 
formula. Choosing for the path (— oo, —e] U eexp([z7r, 0]) U [e, oo) and letting e — > 0, 
we find 



Iiz;r) 



r V 12 ) ./o 



sinh 2xz 



^sinhx sinhrx 
Since the second term is an even function of z, this gives 

- 12^2 + 1 



2z \ dx 

TX / X 



(3.30) 



^{z;t) ^{-z;t) = e 
From (|33C| ) we also get 

/(z + l/2;r) - /(z-1/2;t) = 



24r 



) 



Airiz f°° ( cosh 2xz 



1 \ dx 

TX I X 



(3.31) 



(3.32) 



The integral equals — log(2cos(7rz/T)) (proof left as an exercise). Dividing by 4 and expo- 
nentiating we get the first of the two functional equations 



^(^-l/2;r) 
$(z + l/2;r) 



1 + &{zIt) 



^{z - t/2-t) 
$(z + t/2;t) 



1 + eiz) , 



(3.33) 



and the second can be proved in the same way or deduced from the first using the obvious 
symmetry property 



$(z;t) = $(z/t;1/t). 



(3.34) 



of the function <I>. (Replace x by x/r in ( 3.21 ) 
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4. The functional equations ( |3.33| ) show that 
$(z;r), which in its initial domain of defini- 
tion clearly has no zeros or poles, extends (for 
fixed r with Re (r) > 0) to a meromorphic 
function of z with simple poles at z £ S(r) 
and simple zeros at z G — H(r), where 



Im(z) 
20 r 



10 



■(r) 



T + 



^>0 + 



c 



(3.35) 

In terms of the normalization ( |3.23| ) , this says 
that ^h{z) has simple poles at z £ H(/i) and 
simple zeroes at z G — H(?i), where 

E{h) = (2Z>o + l)iTT + (2Z>o + l)n. (3.36) 



-20 



-10 



-10 



-20L 



10 



' Re(z) 



20 



This is illustrated in Figure ^. Equation ( ^.221) Figure 5: The complex z-plane, showing poles 
expressing $ in terms of the function (x; q)oo (^'s) ^nd zeroes (O's) of $ri(z) at ^ = fe^/^. 
also follows, because the quotient of its left- and right-hand sides is a doubly periodic 
function of z with no zeros or poles and hence constant, and the constant can only be ±1 
(and can then be checked to be +1 in several ways, e.g., by evaluating numerically at one 
point) because the right-hand side of ( 3.22| ) satisfies the same functional equation (3.31) as 
<I>(z;t) by virtue of the Jacobi triple product formula ( p.29| ) and the well-known modular 
transformation properties of the Jacobi theta function. 

5. From ( |3.30| ) we also find the Taylor expansion of I{z; t) at z = 0, 



/(z;r) = 4^Cfc( 



r z 



k=0 

with coefficients Cfc(T) = r~'^Cfc(l/r) given by 
ni 



Co(r) 
C2(t) 
Cfe(r) 



24 

"27' 
2k-i 



Ckir) 







1 



sinhx sinhrx 
for k > 4: even, 

''^^ dx 



dx , 



for k > 3 odd. 



sinhx sinhrx 

By expanding l/sinh(x) and l/sinh(ra;) as power series in and e~'^^ we can evaluate 
the last of these expressions to get 

nfc + l /"OO r) 

k ^ 



Ckir) 



kl 



E 



m,n>0, odd 



{k>3 odd) 



with H(r) as in ( |3.35| ). Dividing by 4 and exponentiating gives the Weierstrass product 
expansion 



$(z;t) = expf^(r + T~i) + Ci{t) z 



TTIZ 



2t 



n 

seH(r) 



S + Z 



~2z/s 



(3.37) 
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of <I>(z;t). From this expansion, one finds that ^(z;t) extends meromorphically to C x 
(C \ (—00, 0]) with simple poles and simple zeros for z € H(r) and z G — S(t) and no other 
zeros or poles. (This analytic continuation can also be deduced by rotating the path of 
integration in ( |3.21| ), e.g. by replacing by /iR(+)/0r for z sufficiently small.) 

6. The quantum dilogarithm is related via the Jacobi triple product formula to the Jacobi 
theta function, which is a Jacobi form, i.e., it has transformation properties not only with 
respect to the lattice translations z z + 1 and z z + t but also with respect to the 
modular transformations r 1— > r + 1 and r 1— > — 1/t. The function $(z;r) has the lattice 
transformation properties ( |3.33|) and maps to its inverse under {z,t) i-^ (z/t, —1/t), but 
it does not transform in a simple way with respect to r 1-^ r + 1. Nevertheless, it has 
an interesting modularity property of a different kind (cocycle property) which is worth 
mentioning here even though no use of it will be made in the remainder of this paper. 
Write (PI) as 

S(z/t:-1M- ln„<Oodd(l + <!"'M--))"' itImM<0, 

where q = e(T). The function S{z;t) has the transformation properties 

S{z; r) = S{z + 1; r) = (1 + q^/^eiz)) S{z + r; r) = S{z + 1; r + 1) = S{z; r + 2) 
amd from these we deduce by a short calculation the two three-term functional equations 

/ 1 N /-2 T t/2 T s , , 

^z-t) = <^>{z±-,T + l)<P{^^,-^) (3.38) 



of ^. This is highly reminiscent of the fact (cf. ||3J]) that the holomorphic function 

Hr) = /(r)-r-27(-l/r), /(r) = 



En>o ^nq'^ if Im (r) > 0, 
E„<o«ng" iflm(r)<0 



associated to a Maass cusp form n(r) on 5L(2,Z) with spectral parameter s, where 
are the normalized coefficients in the Fourier-Bessel expansion of u, satisfies the Lewis 
functional equation 'iP{t) = ^(r-|- 1) + {t+ l)^^*^(:p^) and extends holomorphically from 
its initial domain of definition C\R to C \ (— oo,0]. 

7. Finally, the quantum dilogarithm functions satisfies various five-term relations, of which 
the classical five-term functional equation of Li2(x) is a limiting case, when the arguments 
are non-commuting variables. The simplest and oldest is the identity 

iY;q)^{X;q)^ = {X;q)^ {-YX;qU {Y;q)^ (3.39) 

for operators X and Y satisfying XY = qYX. From this one deduces the "quantum 
pentagon relation" 

^n{p)'^h{q) = ^n{q)'^h{p + q)^h{p) (3.40) 
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for operators p and q satisfying [q,p\ = 2h. Letting Sij be a copy, acting on the i^^ and 
j^^ factors of V (2)V (2)V, of the S-matrix introduced in ( p.l2 ), we deduce from (3.40) the 
operator identity 

S23S12 = S12S13S23 • (3-41) 
It is this very special property which guarantees that the gluing procedure used in the 



definition of (3.16) is invariant under 2-3 Pachner moves on the underlying triangulations 



and produces a true hyperbolic invariant |26, 27, 3C]. This identity is also related to the 



fact (not used in this paper) that the fifth power of the operator which maps a nicely 
behaved function to the Fourier transform of its product with $^(-2) (suitably normalized) 
is a multiple of the identity [^] . 

3.4 A state integral model for Z^'') (M ; h) 



Now, let us return to the analysis of the integral ( 3.16| ) and compare it with the perturba- 
tive SL(2, C) invariant Z^") (M; h, u). Both invariants compute quantum (i. e. /I - deformed) 
topological invariants of hyperbolic 3-manifolds and, thus, are expected to be closely re- 
lated. However, in order to establish a precise relation, we need to face two problems 
mentioned in the beginning of this section: 



i) the integration contour is not specified in ( p. 161) , and 

ii) the integral ( |3.16D does not depend on the choice of the classical solution a. 
These two problems are related, and can be addressed by studying the integral (|3.16| ) in 
various saddle point approximations. Using the leading term in ( 3.26| ), we can approximate 
it to leading order as 

H{M;h,u) ~ / e:^^(Pi'-'P^-i'") dpi . . . dpjv-i , (3.42) 
with the "potential" 

1 ^ 1 
Vip,u) = - J^eiLi2(-exp(<7i(p,2n)+i7rei)) + -/(p,2n,n = 0). (3.43) 



i=l 



As explained below (|3.26| ), the branches of Li2 must be chosen appropriately to coincide 
with the half- lines of poles and zeroes of the quantum dilogarithms in ( 3.16| ). The leading 



contribution to H{M; h, u) will then come from the highest-lying critical point through 
which a given contour can be deformed. 

It was the observation of [^] that the potential V always has one critical point that 
reproduces the classical Chern-Simons action on the "geometric" branch, that is a critical 
point p(^'^°™)(n) such that 

keoM = exp [|-y(p{g--)(n),n)] (3.44) 

and 

^(goom) ^ y (p{gcom) ^ ^ (3 45) 

This identification follows from the fact that both S'q^'^°™^ (u) and matrix elements ( ^.131) - 



(3.14) in the limit h ^ are related to the complexified volume function i(Vol{M;u) + 
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iCS(M; u)). We want to argue presently that in fact every critical point of V corresponds 
to a classical solution in Chern-Simons theory (that is, to a branch of A{1, m) = 0) in this 
manner, with similar relations 



=e^^(p'"'(")'") (3.46) 



and 



5i")(n)=y(p(°)(^x),n) (3.47) 
for some Oi. In particulcir, Iq, (u) as given by (|3.46| ) and m = e" obey (|2.5lD . 



To analyze generic critical points of V, observe that the critical point equations take 
the form 

d ^ d 

2— y(p, u) = -Y, ^^GJi log(l+exp(gi(p, 2u)+i7re^)) + — f{p, 2u, 0) = V j , (3.48) 



where Gji = ■^gi{p,2u) are some constants and the functions ^/(p, 2u,0) are linear. 
Again, the cuts of the logarithm must match the singularities of the quantum dilogarithm. 
Exponentiating ( [3.48 ), we obtain another set of conditions 



rj(x,m) = l Vj, (3.49) 

where 

rJ^ rr^.) = Pvn f 2 V{p,u)) = JJ(1 - e^')'^''^'' CXp ( — 

opj ) ^(ypj 



are all rational functions of the variables xj = e^^ and m = e". Note that an entire family 
of points {p + 27rm | n G Z^~^} maps to a single x, and that all branch cut ambiguities 
disappear in the simpler equations ( |3.49| ). Depending on arg(^) and the precise form of /, 
solutions to ( |3.49D either "lift" uniquely to critical points of the potential V, or they lift to 
a family of critical points at which V differs only by integer multiples of liriu. 

Now, the system ( 0.49 ) is algebraic, so its set of solutions defines a complex affine 



variety 

n = {{xi,. . . ,XN-i,m) £ \rj{xi, . . . ,XN-i,m) = Vj} , (3.51) 

which is closely related to the representation variety C given by A{l,m) = 0. Both gener- 
ically have complex dimension one. Noting that s{xi, . . . ,XN-i,m) = exp (^^) is also a 
rational function, we can define a rational map (j) : by 

(/)(xi, . . .,XN~.i,m) = {s{xi, . . . ,XN-i,m),m) . (3.52) 

The claim in p7| that one critical point of V always corresponds to the geometric branch 



of C means that (/'(T^) (taking an algebraic closure) always intersects C nontrivially, along 
a subvariety of dimension 1. Thus, some irreducible component of (j)ijV)^ coming from 
an irreducible component of IZ, must coincide with the entire irreducible component of C 
containing the geometric branch. Every solution x = x(m) in this component of TZ corre- 
sponds to a branch of the A-polynomial. Moreover, if such a solution x^") (corresponding 
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to branch a) can be lifted to a real critical point p^"^ (u) of V, then one must have relations 
( |3l6D and 



This simple algebraic analysis shows that some solutions of ( 3.49| ) will cover an entire 



irreducible component of the curve C defined by A(l,m) = 0. We cannot push the general 
argument further without knowing more about the reducibility of TZ. However, we can look 
at some actual examples. Computing V{p, u) for thirteen hyperbolic manifolds with a sin- 
gle torus boundary/'' we found in every case that solutions of (3.49) completely covered all 



non-abelian branches a / abel; in other words, (p{TZ) = C , with C = {A{l,m)/{1 — 1) = 0}. 
For six of these manifolds, we found unique critical points p{u) corresponding to every non- 
abelian branch of C at arg(fi) = in. Motivated by these examples, it is natural to state the 
following conjecture: 

Conjecture 2: Every critical point of V corresponds to some branch a, and all a ^ abel 
are obtained in this way. Moreover, for every critical point p*-") (u) ( corresponding to some 
branch a) we have ( 3.46| )- (|3.47| ) and to all orders in perturbation theory: 



N N-l 



z("\M;h,u) = V2 [ Y[<^n{9^{p,2u)+e,ii7^ + h)y'e^J(P''-'^)- 1]^%' (3-53) 

where Ca is an arbitrary contour with fixed endpoints which passes through p(") (u) and no 
other critical point. 

A slightly more conservative version of this conjecture might state that only those a 
that belong to the same irreducible component of C as the geometric branch, a = geom, 
are covered by critical points of V. Indeed, the "abelian" branch with /abel = 1 is not 
covered by the critical points of V and it belongs to the separate component (/ — 1) of the 
curve A{l,m) = 0. It would be interesting to study the relation between critical points 
of V and irreducible components of C further, in particular by looking at examples with 
reducible A-polynomials aside from the universal (/ — 1) factor. 

The right-hand side of (|3.53D is the proposed state integral model for the exact per- 
turbative partition function of 5L(2,C) Chern-Simons theory on a hyperbolic 3-manifold 
M with a single torus boundary S = T^. (A generalization to 3- manifolds with an arbi- 
trary number of boundary components is straightforward.) This state integral model is a 
modified version of Hikami's invariant (3.16). Just like its predecessor, eq. ( p. 53 ) is based 



on an ideal triangulation {Aj}^^ of a hyperbolic 3-manifold M and inherits topological 
invariance from the pentagon identity ( 3.40]) of the quantum dilogarithm. 



However, in writing ( |3.53| ) we made two important modifications to Hikami's invariant 



(13.161) . First, we introduced contours Cq running across the sadle points Pa, which now 



encode the choice of a classical solution in Chern-Simons theory. Second, in ( |3.53| ) we 



^'^Namely, the complements of hyperbolic knots 4i(k2i), 52(k32), 12n242 (3i, (-2,3,7)-pretzel knot), 
6i(k4i), 63(k643), 72(k42), 73(k52o), 74(k628), 10i32(K59), 10i39(K522), and lln38(K5i3), as well 
as the one-punctured torus bundles L^R and LR"^ over (also knot complements, but in a manifold other 
than S^). 
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introduced an extra factor of VSirhe"'^, which is needed to reproduce the correct asymptotic 
behavior of Z^°'^ (M; h, u). To understand this correction factor, we must look at the higher- 
order terms in the expansion of Z(")(M;/i,u). By using ( |3.26[ ) , one can continue the saddle 



point approximations described above to arbitrary order in h. The result has the expected 
form (|r|), 

Z(")(M; h, u) = exp {\s^^\u) - ^-J^'^ log^ + S^^'liuW] , (3.54) 

\ n=0 / 

with the correct leading term Sq"^ (n) that we already analyzed above, cf. eq. ( p.47| ) . 

Let us examine the next-leading logarithmic term. Its coefficient (5^") receives contribu- 
tions from two places: from the prefactor {4TTh)~^^~^^^'^ in ( 3.53| ), and from the standard 



Gaussian determinant. The former depends on the total number of tetrahedra, N, in 
the triangulation of M and therefore must be cancelled (at least partially) since the total 
integral ( |3.53| ) is a topological invariant and cannot depend on A^. This is indeed what 
happens. For example, in a saddle point approximation around a nondegenerate critical 
point p^"\u), the contribution of the Gaussian determinant goes like ~ and ex- 

actly cancels the contribution of the prefactor ~ h~^^~^^^'^ . An example of such critical 
point is the critical point p^^'^°^\u) corresponding to the geometric branch. Therefore, 
the asymptotic expansion of the integral ( 3.53| ) around the critical point pfe^°™)(u) has the 



form ( 3.54| ) with 

^(geom) ^ Q ^ (3 55) 

which is the expected result. Indeed, as explained e.g. in ||2^, the rigidity of the flat 
connection associated with a hyperbolic structure on M implies = = 0, so 

that (13) gives jCs^"'") = 0. 

Using ( |3.26 ) , one can also calculate the higher-order perturbative coefficients sii'\u), 



with n > 1. In the following section, we carry out this analysis to high order for the figure- 
8 knot complement and find perfect agreement with the results obtained by methods of 
Section ^. (Other interesting examples and further checks will appear elsewhere p^.) 

Note that the S^\u)^s do not depend on the details of the contours Ca- The only 
part of ( |3.53| ) which actually depends on the details of Ca is exponentially suppressed and 
is not part of the perturbative series ( |3.54 ). Finally, we also note that we use only those 



critical points of the full integrand ( p. 53 ) which correspond to critical points of V in the 
limit h ^ 0. For any fixed ^ > 0, the actual integrand has many other critical points 
which become trapped in the half-line singularities of quantum dilogarithms as — > 0, so 
the integrals over them do not have a well-behaved limit. 



^^Recall that throughout this paper Z^''\M; K) (resp. Z'"'(M; h,u)) stands for the unnormalized pertur- 
bative Gc invariant. A normalized version, obtained by dividing by ^(S^), has an asymptotic expansion of 
the same form (resp. (|^)) with the value of 5'"' shifted by dim(G) for every a. This is easy to see 



from (2.18) 
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We conclude this section by observing that Conjecture 2 implies Conjecture 1. From 
(3.26), we can write an asymptotic double series expansion (for very small p): 



Mpo +p) = exp (f^Bnil + 4) Li2-n(-e^'») 



^2 2h^ ^ ^ n! 

\n=0 



oo oo 



exp 



EE °y|'ff,f ".-.->(-^°)ftV|. (3.56) 



Using this formula and taking into account the shifts by ib(i7r + h), we expand every 
quantum dilogarithm appearing in the integrand of ( 3.53| ) around a critical point p^"^. At 



each order in h, the state integral model then reduces to an integral of a polynomial in p 
with a Gaussian weight. Due to the fact that Li^ is a rational function for A; < 0, the 
coefficients of these polynomials are all rational functions of the variables x^") = exp(p(")) 
and m. Therefore, the resulting coefficients sii'\m), for n > 1, will also be rational 



functions of x^") and m. At m = 1, the solutions x(")(m) to the rational equations ( 3.49] ) 



all belong to some algebraic number field IC C Q, leading immediately to Conjecture 1. 
In particular, for the geometric branch a = geom, the field IK is nothing but the trace 
field Q(trr). 



4. Examples and computations 

In the previous sections, we described several ways to approach perturbative Chern-Simons 
theory with complex gauge group. Our discussion was guided by the aim of computing the 
all-loop perturbative partition function Z^^\M;h) in the background of a nontrivial flat 
connection A^^'' labeled by a homomorphism p : tti{M) Gc- 

In this section, our goal is to illustrate these methods in concrete examples by cal- 
culating perturbative invariants Sn^ (also denoted S^\u)) to high loop order. Which 
3-manifolds shall we choose for our examples? On one hand, the approach based on quan- 
tization and analytic continuation most directly applies to 3-manifolds with boundary, such 
as knot complements. (The invariants of closed 3-manifolds without boundary can be ob- 
tained by gluing two manifolds, M+ and M_, along a common boundary S, as in ( |2.19| ).) 
On the other hand, the state integral model described in Section ^ does not require M 
to have a boundary, but in the present form applies only to hyperbolic 3-manifolds. (We 
believe this state integral model can be extended to arbitrary 3-manifolds, but we do not 
pursue it here.) Therefore, if we wish to compare the results of different methods, it is 
convenient to choose examples of 3-manifolds that are knot complements and hyperbolic 
at the same time. In other words, we take M to be a complement of a hyperbolic knot K 
in the 3-sphere: 

M = S^\K. (4.1) 

Among hyperbolic knots with a small number of crossings, the simplest ones are the figure- 
eight knot (also denoted 4i) and the three- twist knot (also denoted 62) shown in Figure |6|. 
Among all hyperbolic knots, the complement of the figure-eight knot has the least possible 
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volume Vol(4i) = 2.0298832128..., whereas the volume of the 62 knot complement is 
Vol (52) = 2.8281220883.... In both cases, the knot group 7ri(M) is generated by two 
elements a and b, such that a~^bab~^a = ba^^bab^^ for the figure-eight knot and waw~^ = 
b, with w = a^^ba^^b^^ab^^ , for the 52 knot. For the figure-eight knot, the trace field 
Q(trr) is the imaginary quadratic field Q(\/— 3), while for the 62 knot it is a cubic field 
with one complex place and discriminant —23. 




b) 




Figure 6: a) figure-8 knot 4i, and b) three-twist knot 52- 



we fix the homomor- 



This section consists of two parts. In the first part (Section |4.l| 
phism p = geom and compute the perturbative SL{2, 
hyperbolic knot complements, M = \ 4i and M = \ 52- In particular, we shall verify 



invariants si^^"™^ for the simplest 



the arithmeticity conjecture formulated in Section 2.1. In the second part (Section [4. 2]) , we 
look more closely at the perturbative invariants of the figure-eight knot complement and 
study how the Sn^^s depend on the homomorphism p (equivalently, we study each sli^\u) 
as a function of u and a). 



4.1 Numerical verification of the arithmeticity conjecture 

In this section, we explicitly calculate the perturbative SL{2,C) invariants S^\u) for the 
simplest hyperbolic knot complements and for the special value of n = (corresponding 
to the discrete faithful representation p = geom). It was conjectured in |2^ that such 
invariants (in fact, for arbitrary value of u) can be extracted from the asymptotic expansion 
of knot invariants computed by Chern-Simons theory with compact gauge group G in the 
double-scaling limit ( |1.14 ). In the present case, G = SU{2) and the relevant knot invariant 
is 



where Jn{K] 



^^^'^^ ' Jjv (unknot; e2W^) ' 
is the 'W-colored Jones polynomial" of normalized such that 



Jjv (unknot; q) 



-rN/2 



-N/2 



7I/2 



-1/2 



(4.2) 



(4.3) 



is the "quantum dimension" of the A^-dimensional representation of SU{2). As explained in 



the Introduction and then in more detail in Section 2.2.1 , the A^-colored Jones polynomial 
J]\f{K;q) is computed by the Feynman path integral in Chern-Simons gauge theory on 
with a Wilson loop operator, Wji{K), supported on K and labeled by the A^-dimensional 
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representation R of SU{2), cf. ( 1.13 ). According to the Generalized Volume Conjecture p2| 
the invariant J^{K; N) should have the asymptotics^^ 



as ^ 00, where 



n=0 



/27rA 



n— 1^ 



(4.4) 



(4.5) 



f(Si X D^; ivr/A^) 

5fe--)(0) = s„-2"-i (n/1) 
5fe--)(0) = si + ilog2 (n = l) 

are the perturbative 5L(2,C) invariants of M = S'^ \ AT. (Here, we view the solid torus, 
X D'^, as the complement of the unknot in the 3-sphere.) Specifically, we have 



so = i(Vol(M) + iCS(M)) 



(4.6) 



(Volume Conjecture) and si is the Ray-Singer torsion of M twisted by a flat connection, 
cf. eq. (|2.3|). The arithmetic TQFT conjecture described in Section [2?^ predicts that 



si G 



logK, 



s„ G IK (n > 2), 



(4.7) 



where K is the trace field of the knot. In this subsection, we will give numerical computa- 
tions supporting this conjecture for the two simplest hyperbolic knots 4i and 52- 

The formulas for J{K] N) in both cases are known explicitly, see e.g. We have 

N-l 

J(4i;iV) = ^(g)™((?-i)^, (4.8) 
m=0 

N-l m 

J{W,N) = Y.1~^"''''^'(^^'rn/iQ~')k, (4.9) 

m=0 k=0 

where q = e^'^*/^ and {q)m = (1 ~ f?) " ' ' (1 ~ q"^) is the t^-Pochhammer symbol as in 
Section The first few values of these invariants are 



N 


:^(4i;A) 




1 


1 


1 


2 


5 


7 


3 


13 


18 -5q 


4 


27 


40 - 23q 


5 


44 - Aq"^ - Aq^ 


46 - 55q - 31g2 + 


6 


89 


120 - 187g 


7 


100 - 14g2 - 25q^ - 25g^ - Uq^ 


-154^ - 88^2 _^ 47^3 _^ 58^4 _^ 77^5 


8 


187 - i5q - 45^3 


-84 - 407g - 150g2 + 96^^ 



'Note that in our notations e^" = q = e 



2-iTi/N 
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Using formula ( |4.8| ) for N of the order of 5000 and the numerical interpolation method 
1], we computed the values of s„ for < n < 27 to very high 



explained in |56] and || 
precision, finding that ( [4. 4]) holds with 



So 



27r2- 



si 



logs, 



(4.10) 



the first in accordance with the volume conjecture and the second in accordance with the 
first statement of ( [4.7D , since IC = Q(\/— 3) in this case, and that the numbers 



s'^ = sn- (6^^)" ' (n > 2) (4.11) 
were very close to rational numbers with relatively small and highly factored denominators: 



n 


2 


3 


4 


5 


6 


7 


8 


9 


S'n 


11 

12 


2 


1081 
90 


98 


110011 
105 


207892 
15 


32729683 
150 


139418294 
35 



n 


10 




27 


s' 


860118209659 




240605060980290369529478710291172763261781986098552 


10395 


814172781296875 



confirming the second prediction in ( |4.7| ). Here D{z) is the Bloch-Wigner dilogarithm (p.4|). 
(That Sn is a rational multiple of (\/— 3)" ^, and not merely an element of Q(\/— 3) is a 
consequence of (|4.4| ) and the fact that ^7(41; A^) is real.) 

Actually, in this case one can prove the correctness of the expansion rigorously: the 
two formulas in ( [4.1C| ) were proved in Q and the rationality of the numbers s'^ defined 
by (|13) in |lO|, |l|] and [||], and can therefore check that the numerically determined 
values are the true ones (see also [^] for a generalization of this analysis). In the case of 
52, such an analysis has not been done, and the numerical interpolation method is therefore 
needed. If we try to do this directly using eq. ([4.9|), the process is very time-consuming 
because, unlike the figure-8 case, there are now 0{N'^) terms. To get around this difficulty, 
we use the formula 



E 



iQ)k 



iQ)r 



q 



fc=0 ^^'^ 



(4.12) 



k=0 '•'^"^ k=0 '-'^'k 

which is proved by observing that both sides vanish for m = — 1 and satisfy the recursion 
= (1 — (f^)tm-i + (f^ /{Q)m ■ This proof gives a way to successively compute each tm 
in 0(1) steps (compute y = q"^ as q times the previous y, {q)m as I — y times its previous 
value, and then tm by the recursion) and hence to compute the whole sum in ( |4.9| ) in only 
O(A^) steps. A short PARI program for this is 

{f(x) = q=exp(2*Pi*I*x) ;y=l;p=l;s=l;t=l; 
f or (n=l ,denom(x) -1 ,y*=q;p*=l-y ; t+=y"n/p"2 ; s+=p"2*conj (p*t) ) ; s} 



and is very much faster than the computation based on the original formula. For instance, 
to compute J {52] N) to 200 digits for N = 100 and N = 1000 using (U) takes 1.6 seconds 
and 262 seconds, respectively, whereas with the improved algorithm these times are reduced 
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to 0.03 seconds and 0.4 seconds. The interpolation method can therefore be carried out to 
just as high precision as in the figure-8 case. 

The results are as follows. The first coefficient is given to high precision by 

so = (Li2(a) + ^ log(a) log(l - a)) + | , (4.13) 

in accordance with the prediction (4.6), where a = 0.87743 • • • — 0.74486 . . . i is the root 
of 

- + 1 = (4.14) 

with negative imaginary part. The next four values are (again numerically to very high 
precision) 

1 1 + 3a 198a2 + 1452a - 1999 

'^ = 4^°^^^' = 24T23^ ' 

_ 465a2 - 465a + 54 _ -21033020^ + 55115a + 5481271 

~ 2^233 ' ~ 240 • 235 ' 

in accordance with the arithmeticity conjecture since K = Q(a) in this case. 

These coefficients are already quite complicated, and the next values even more so. 
We can simplify them by making the rescaling 

4 = 5nA"-i (4.15) 

{i.e., by expanding in powers of 2'iri/XN instead of Ini/N), where 

X = a^{3a-2f = a-^{a'^ -3f . (4.16) 

(This number is a generator of p^, where p = (3a — 2) = (a^ — 3) is the unique ramified 
prime ideal of K, of norm 23.) We then find 

4 = -j3 (12 - 19 a + 86), 

4 = -| (2a2 + 5a -4), 

s'^ = ^ (494 a^ + 12431 a + 1926) , 

s'g = - i (577 a^ - 842 a + 1497) , 

4 = Ygigo (176530333 a^ - 80229954 a - 18058879), 

s'j = -220 (99281740 a^ + 40494555 a + 63284429) , 

4 = -io^bo (3270153377244 a^ - 4926985303821a - 8792961648103), 

s'g = (9875382391800 a^ - 939631794912 a - 7973863388897) , 

•^'lo = - 15966720 (188477928956464660 a^ + 213430022592301436 a + 61086306651454303) , 
•s'n = - 1209600 (517421716298434577 a^ - 286061854126193276 a - 701171308042539352), 

with much simpler coefficients than before, and with each denominator dividing (n + 2)! . 
These highly nontrivial numbers give a strong experimental confirmation of the conjecture. 

We observe that in both examples treated the first statement of the conjecture ( [4.7| ) 
can be strengthened to 

exp(4si) G K. (4.17) 
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It would be interesting to know if the same statement holds for all hyperbolic knot comple- 
ments (or even all hyperbolic 3-manifolds). Another comment in this vein is that ( |4.6| ) also 
has an arithmetic content: one knows that the right-hand side of this equation is in the 
image under the extended regulator map of an element in the Bloch group (or, equivalently, 
the third algebraic K-group) of the number field K. 



4.2 Perturbative SL{2, C) invariants for the figure-8 knot complement 

In this subsection, we explicitly calculate the perturbative SL{2, C) invariants Sn^ for 
the figure-8 knot complement, first via quantization of the moduli space of flat SL{2, C) 
connections and then, independently, using the state integral model formulated in Section |3|. 
We find perfect agreement between the results, thus, supporting our proposal ( 3.53| ). 



Let K be the figure-eight knot (shown in Figure 0.a) and let M be its complement 
in the 3-sphere. In order to compute the perturbative invariants Sn for an arbitrary 
p : 7ri(M) SL(2, C) we first need to review the classical geometry of M in more detail 
and, in particular, to describe the moduli space of flat SL{2, C) connections on M. As 
we already mentioned earlier, the knot group 7ri(M) is generated by two elements, a and 
b, such that a~^bab~^a = ba~^bab~^ . The corresponding representation into SL{2,C) is 
given by 



p{a) =i 1 , p{b) =i 1 , (4.18) 





where C = (—1 -|- \^—3)/2 is the cube root of unity, = 1- 

The complement of the figure-eight knot can be also represented as a quotient space 
IHI'^/r (p.6|), where the holonomy group F is generated by the above two matrices. Specifi- 
cally, we have 

T^PSL{2,Ok), (4.19) 

where Ok is the ring of integers in the imaginary quadratic field K = Q{\/ —3). The 
fundamental domain, J^, for F is described by a geodesic pyramid in with one vertex at 
infinity and the other four vertices at the points: 



Explicitly, we have 



vi= J 



Vo = — 

2 



73. 

6 ^ 



3^ 



(4.20) 



^^3 



1 \/3. 

- + —I + 

2 6 



1 



^/3 



i + 



J" = {z + G I z G jr^, xl + |z|^ > 1} , 



(4.21) 
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where z = xi + 1x2 and 

jr = {z G C I < Re(z), ^Re (z) < Im(z), Im (z) < ^(1 - Re(z))} 

v 3 \/3 

U{z e C I < Re(z) < J, - ^Re(z) < Im(z) < ^Re(z)} (4.22) 

2 y 3 y 3 

is the fundamental domain of a 2-torus with modular parameter t = Q. The region of 
large values of x^ in T corresponds to the region near the cusp of the figure-eight knot 
complement M. 

The standard triangulation of the figure eight knot complement comprises two ideal 
tetrahedra of opposite simplicial orientations, as in Figure glued together in the only 
nontrivial consistent manner possible, 

M = U . (4.23) 

Here, z and w are complex numbers, representing the shapes of the ideal tetrahedra; we 
take A^ to be positively oriented and A^ to be negatively oriented. As explained in 
Section 3.1, the shape parameters z and w must obey consistency relations, which in the 
case of the figure-eight knot reduce to a single algebraic relation (see e.g. Chapter 4 of |46| 
and Section 15 of @): 

{z - l){w - 1) = z'^w'^ . (4.24) 

The shape parameters z and w are related to the SL(2,C) holonomy eigenvalue, I, along 
the longitude^'' of the knot in the following way 

^ = -/, -!^ = -l-\ (4.25) 
z — 1 w — 1 



which automatically solves the consistency condition (4.24). Similarly, the holonomy eigen- 
value m = e" around the noncontractible meridian of the torus is given by 



zw = m. (4.26) 

Given the decomposition (|4.23| ) of the 3-manifold M into two tetrahedra, we can find 
its volume by adding the (signed) volumes of A^ and A^;, 

Vol (M; u) = Vol (A^) - Vol (A^) = D{z) - D{w) , (4.27) 

where the volume of an ideal tetrahedron is given by (^^). Similarly, following the pre- 



scription in Section 3.1, the complexified volume can be given by^^ 

i(Vol (M; u) + iCS(M; u)) = L{z; 0, 0) - L{w; 0, 0) - vu - mu . (4.28) 

Notice that due to the consistency relation (4.24) the total volumes ( [4.27| ), ( 4.28]) are 
functions of one complex parameter, or, equivalently, a point on the zero locus of the 
A-polynomial, A{l,m) = A{—e^,e^) = 0. 



1-cycle onT, — which is contractible in the knot complement M. 
^^This expression differs slightly from the one given in |40|, because we use 2v + 2Tvi rather than 2v as 



the logarithm of the "longitudinal" holonomy, as mentioned in Footnote hsi of Section 



3.1 
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The classical A-polynomial of the figure-eight knot is 

A{1, m) = {l- l)(m^ - (1 - - 2m^ - + m^)l + rriH"^) . (4.29) 

Observe that A{1~^ ,m~^) ~ A{l,m) and that A{l^^,m) ~ A{l,m), the latter reflecting the 
fact that the figure-eight knot is amphicheiral. The zero locus of this A-polynomial has 
three branches: the abelian branch ^abei = 1 (or t;abei = — ^vr) from the first factor, and two 
other branches from the second, given explicitly by 

1 — m? — 2rn^ — + rrfi 1 — 
/geom,conj (m) = ± A(m) (4.30) 

or 

'l'geom,conj (^) — log( Zgeom,conj (c )) i (^-^l) 

where we have defined 



A(m) = z\/-m-4 + 2m-2 + 1 2m? - . (4.32) 

Note that A(m) is an analytic function of m = around u = 0, as are l{m) and v{u). In 
particular, A(l) = \/— 3 is the generator of the trace field IK = Q(trr) for the figure-eight 
knot. 



Since there are only two non-abelian branches, they must (as indicated in Section 2.3) 
be the geometric and conjugate ones. Alternatively, eliminating shape parameters from 
the (geometric) eqs.(4.26) and (4.25) implies precisely that the non-abelian factor of the 
A-polynomial vanishes. As expected for an amphicheiral knot, due to the symmetries 
discussed in Section (2.3), Vgeom{u) = —Vconiiu)- The exact labelling of branches (geometric 
vs. conjugate) can be verified by computing the hyperbolic volume Vol (M; u) of the figure- 
eight knot complement via an ideal triangulation, as in ( 4.28| ), and comparing the result 
to the integral ( |2.52 ) on the geometric branch. One obtains precisely 

S^^^^'^^u) = ^(Vol (M; 0) + iCS(M; 0)) + P {vgeom{u') + i7r)du' 
^ Jo 

= -(Li2(z) - Li2(w) + -log2:log(l - z) - -log(w) log(l - w) + uv + ittu) 

= ^(Vol (M; u) + iCS(M; u)) + Vgeom{u)Re {u) + iiru , (4.33) 

confirming ( ^153] ). Note that CS{M;u) = for any imaginary u. 

From the relation ( 4.25| ) we find that the point {l,m) = (—1, 1) corresponding to the 
complete hyperbolic structure on M is characterized by the values of z and w which solve 
the equation 

z"^ -z + l = 0. (4.34) 

In order to obtain tetrahedra of positive (signed) volume, we must choose z to be the root 
of this equation with a positive imaginary part, and w its inverse: 

z = , w = . (4.35) 

These value correspond to regular ideal tetrahedra, and maximize (respectively, minimize) 
the Bloch-Wigner dilogarithm function D{z). 
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4.2.1 Quantization 

In order to test the proposed state integral model ( p. 531 ), we first compute the perturbative 
SL{2,C) invariants of M via quantization of A^fiat(G'C) 5]), as outlinted in Section |2.2| . 
Recall the expression (1^) for the perturbative Chern-Simons partition function, 

Z(")(M; h, u) = exp [\s'^^\u) - ^<5(-) log^ + f; S^:i^{u)tA , (4.36) 

\ n=0 / 

where for the figure-eight knot complement a G {geom, conj, abel}, and, as explained e.g. 
in H, H], ^("T^abei) ^ 

In Section [2^ we saw how quantization of the relation A{1, m) = leads to the equation 

A{i, m) (M; ^i, u) = , (4.37) 

which then determines the perturbative Gc invariant Z^'^\M;h,u). As explained earlier, 
cf. eqs. ( |1.17 ) and ( ^.36 ), a simple way to find the operator A{l,7n) is to note that it also 
annihilates the polynomial invariants (of the knot K) computed by Chern-Simons theory 
with compact gauge group G. In the present case, G = SU{2) and the equation ( 2.36| ) 
takes the form of a recursion relation on the set of colored Jones polynomials, {Jn{K; q)}. 
Specifically, using the fact that / acts by shifting the value of the highest weight of the 
representation and writing A{l,7fi) = Yl'j=Q'^j{''^^QW ^ ™ ( 2.43| ), we obtain 

d 

^a,((?"/2,g)J„+,(K;g) = 0. (4.38) 
i=o 

It is easy to verify that the colored Jones polynomials of the figure-eight knot indeed satisfy 
such a recursion relation, with the coefficients (see also^^ |17, p!8|): 

2 

ao{m, q) = — :rj7 , (4.39a) 

, , , 1 + (g^ - 2q)rh? — (q^ — q"^ + q)m'^ — {2q^ — q^)m^ + q^rn^ . ^ , , 

ai [m, q) = 1/2-2/1 ^2-2 ^4 3^^6^ ' (4.39b) 

1 - (2g2 - q)m^ - (q^ - + q^)m^ + (g^ - 2q^)m^ + q^m^ 
02 m, q) = r^— ^3 5-74 g^rgT , 4.39c) 

q TTi 

asim, q) = 1/2/1 ^ 2 - 2\f 1 , 5 - 4\ ' (4.39d) 

g^/^(l + 1 + q°m^) 

Note that in the classical limit q = e^^ — > 1, we have 

I(/,m) ^. (4.40) 

^' ' h-^o m2(m2 - l)(m2 + 1)2 ^ ^ 



^■^Note that |l^] look at asymptotics of the colored Jones polynomial normalized by its value at the 
unknot, while the SL{2, C) Chern-Simons partition function should agree with the unnormalized colored 
Jones polynomial. Hence, we must divide the expressions for aj there by (m^g-'''^ — 9 "'''■^) to account for 
the difference, introducing a few factors of q^^^ in our formulas. 
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We could multiply all the a^-'s by the denominator of ( 4.40[) (or a g-deformation thereof) to 



obtain a more direct correspondence between A{l,7h) and A{l,m), but this does not affect 
any of the following calculations. 

To find the exact perturbative invariant Z^°'\M;h,u) for each a, we reduce rh to 
a classical variable m (since in the u-space representation the operator rh just acts via 
ordinary multiplication), expand each of the above Oj's as aj{m,q) = X^^o ^i>p("'')^^' ^^'^ 
substitute the aj^p into the hierarchy of differential equations derived in Section |2.2| and 
displayed in Table |l|. The equations are then solved recursively on each branch a of the 
A-polynomial to determine the coefficients sii'\u). 

Geometric branch 

For the geometric branch, the solution to the first equation in Table ^ is chosen to be 

5fe<^°-)(u) = Vol (4i; 0) + / duv^eomiu) + iTTU , (4.41) 

^ Jo 



with Vgeom{u) as in ( [4.31 ). The integration constant z(Vol(0) + iCS(O)) is not important 



for determining the remaining coefficients (since only derivatives of Sq appear in the equa- 
tions), but we have fixed it by requiring that 4^°°"'Ho) = |(Vol(4i;0) +iCS(4i;0)) = 



I Vol (4i; 0) = (1.01494 . . as expected for the classical action of Chern-Simons theory. 

Substituting the above Sq^^°^\u) into the hierarchy of equations, the rest are readily 
solved^^ for the subleading coefficients. The first eight functions S^^°^^ (u) appear below: 



Si{u) = --log 



2 ° V 2 

S2{u) = — ^ ^ il-v? - 2m'^ + Ibm^ - 2m* - + m^^) , 

12A(m)'^m° 

S'iiu) = ^ ^ (l-m^ - 2m^ + 5m^ - 2m^ - m^^ + m^^) , 

SAu) = — , - 4m2 - 128m^ + 36m^ 

+1074m* - 5630m^° + 5782mi2 _^ 7434^14 _ issilm^^ + 7484m^^ 

+5782m20 - 5630m22 + I074m'^'^ + 36m'^^ - 128m'^^ - Am^^ + m^'^) 
2 

Sbiu) = (1 + bm^ - 35m^ + 240m^ - 282m* - 978m^° 

3A(m)^^m^* 

+3914mi2 - 3496ml'' - 4205mi6 + 9819mi* - 4205m20 - 3496m22 
+3914m24 - 978m2'5 - 2%2m^^ + 240m3° - 35m32 _^ 5^34 _^ ^36^ ^ 

'56(^) = n.^.Tus 28 + 2"^' + 169m^ + 4834m6 
945A(mj^°m^^ 



■^^It is computationally advantageous to express everything in terms of m = e" and rn-^ = etc. 
when implementing this on a computer. 
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-24460m^ + 24U72m^° - 65355^^2 _ 3040056mi^ + 13729993771^*^ 
-15693080mi^ - 36091774^2° + 129092600m22 - 103336363m24 
-11971571677726 + 27078556577728 - 119715716777^0 - 103336363777^2 
+12909260077734 - 36091774777^6 - 15693080777^^ + 13729993777^° - 3040056777 
-65355777^4 + 241472777^6 - 24460777^8 + 4834777^° + 169777^2 _^ 2777^^ _^ ^56^ ^ 

'^7(") = ,^^^_^u8^3o (^ + ^''"^^ ~ l'^^"^^ + 3373777^ + 9683777^ 



45A(777) 

-II6636777IO + 562249777^2 - 515145777^4 - 3761442777^6 _^ 1493987I777I8 
-1552311777720 - 2906145877722 + 9645533577724 - 7152226177726 
-8092952277728 + 17907431577730 - 80929522777^2 - 71522261777^4 
+96455335777^6 - 29061458777^8 - 1552311777740 + 1493987177742 
-376144277744 _ 515145^46 _^ 5622497^48 - 116636777^0 

+9683777^2 _^ 3373^54 _ 175^56 _^ 47^58 _^ ^60^ ^ 

Ssiu) = , / ( 1 + 447772 _ ggg^4 

^ ^ 9450A(777 ) 2 l77740 ^ 

-257567776 + 253397778 - 2848194777^0 - 28212360777^2 
+216407820777^4 _ 1122018175771^6 _ 266877530777 
+1913404485277720 - 7657153250277722 + 75899475728777 
+32445443882877726 - 120620690118277728 + 1153211096310777 
+1903970421177777^2 - 5957756639958777^4 _^ 4180507070492777^6 
+4649717451712777^8 - 101323727219497774° + 4649717451712777 
+418050707049277744 - 595775663995877746 + 1903970421177777 
+1153211096310777^0 - 1206206901182777^2 _^ 324454438828777^4 
+75899475728777^6 _ 76571532502777^8 _^ 191340448527776O 
-26687753O77762 - II22OI8I7577764 + 21640782077766 - 2821236O77768 
-2848194777™ + 25339777^2 _ 25756777^^4 _ QgQ^7e + 44,^78 _^ ^80^ ^ 

Table 2: Perturbative invariants S^°"^\u) up to eight loops. 



as 

,24 

,30 



42 

48 



According to (p.3|), the coefficient S^^ ^u) in the perturbative Chern-Simons par- 
tition function should be related to the Reidemeister-Ray-Singer torsion of M twisted by 
^(geom)^ which has been independently computed. Our function matches24 that appearing 

■^^To compare with jisj, note that fcthorc = fchero ~ ivr/ft, and uthorc ~ 2uhGrc- The shift by — log(7r) is 
directly related to a jump in the asymptotics of the colored Jones polynomial at it = 0. 
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in e.g. |23], up to a shift by — logvr. The constants of integration for the remaining coeffi- 
cients have been fixed by comparison to the asymptotics of the colored Jones polynomial, 



using ( |4.5| ) and the results of Section |4J . 
Conjugate branch 

For the "conjugate" branch, the solution for Sq{u) is now chosen to be 



-Vol (4i;0) + 



du V, 



conj 



[u) + inu (mod 2-ku), 



(4.42) 



so that S, 



(conj) 
' 



-Sq^^"^\u). As for the geometric branch, this is then substituted 



into remainder of the hierarchy of equations. Calculating the subleading coefficients, the 
constants of integration can all be fixed so that 



5(-°J)(n) = i-ir+'sl^ir\^), (4.43) 

This is precisely what one expects for an amphicheiral knot when a conjugate pair of 
branches coincides with a "signed" pair, as discussed in Section 2.3. 

Abelian branch 

For completeness, we can also mention perturbation theory around an abelian flat connec- 
tion ^{^t)*^!) on M, although it has no obvious counterpart in the state integral model. 
For an abelian flat connection A^^^^^\ the classical Chern-Simons action ( [l.g| ) vanishes. 



(4.44) 



This is exactly what one finds from (|2.52| ): 



(abcl) 



(infabei(^) + i'^u = (mod 2Triu) 



fixing the constant of integration so that Sq^^'^^^ (0) = 0. From the hierarchy of differential 
equations, the first few subleading coefficients are 



i(abel) 



5. 



(abcl) 



S. 



(abcl) 



(n 



log 



m[m 



I) 



4(m^ 



(1 — 'im? + m^) 



4(m^ 



(41 - 656m2 + 4427m^ - 16334m^ + 35417m^ - 46266m^° 



3(1 - 3m2 + m4) 
-]-35417mi2 _ i6334m^^ + U27m^^ - 656m^^ + 41m'^°) 



J"^'')(n) =0 



Table 3: Perturbative invariants sii'^^^\u) up to six loops. 
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In the language of Section |2.3| , the abelian branch must be its own "signed pair," guaran- 
teeing that all even S2^^^\u) vanish. 

Having determined the perturbative coefficients for the three branches of classical so- 
lutions (in principle to any order desired) , we can now compare them with the computation 
of the state integral model. 

4.2.2 State integral model 

The state integral model ( |3.53| ) for the figure-eight knot complement gives: 



dp 



<i>n{p + iTT + h) 



Jc^ '^h{-p -2u-iTT-h) 



^-'j^u{u+p)~u 



(4.45) 



There are two tetrahedra (A^ = 2) in the standard triangulation of M, and so two quantum 
dilogarithms in the integral. There is a single integration variable p, and we can identify 
giip,u) =p, g2ip,u) = -p-2u, and f{p,2u,h) = -4:u{u+p). 




Figure 7: Plots of \e'^^^'P''^^ and its logarithm at u = and h — 



It will be convenient here to actually change variables p^p — u — in — h, removing the 
{iTT + H) terms in the quantum dilogarithms, and obtaining the somewhat more symmetric 
expression 



1 



V2nh 



e li 



27Tiu _ 

e ft 



dp 



^h{-p - u) 

Co. 



2pu 



(4.46) 
(4.47) 



We define e'^^^'P'") 



*ft(-p-u) 



e-TT . Figures0and|showplotsof |e^(^'P'")| and log |e^(^'P'")| 



Re T{h,p,u) at fi = i/?> and two values of u. The half-lines of poles and zeroes of the two 
quantum dilogarithms combine into similar singularities for T(h,p,u), as is depicted in 
Figure ^; note the splitting of these poles and zeroes by an amount 2u. 

After our change of variables, the "potential" function V{p,u) as in ( p. 43 ) is now seen 
to be 

V{p, u) = l [Li2(-eP-") - Li2(-e-P-") - 'ipu + Airiu] . (4.48) 
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Figure 8: Plots of je^C^'P'")! and its logarithm at u — and h — 



Instead of looking directly at 
equation 

r(x, m) 



-y = to find its critical points, we consider the simpler 



e 9p 



V 



-T7 ^7^ V = 1 ' (4-49) 

m'^ym + x)[\ + mx) 

in terms of x = and m = e". This clearly has two branches of solutions, which both lift 
to true critical points of V, given by 



P 



(geom,conj) 



u) = log 



1 



m 



^ m^A(m) 



2m? 



with A(m) defined in ( 4.32 ). The lift is unique if (say) i 
correspond to the geometric and con- 
jugate branches of the A-polynomial, 
which can be verified by calculating^^ 
m) = exp {-§^V{p, u)) = [ 



(4.50) 

>o. We claim that these 



Im(p) 
15f 



m+x 



1 1/2 



and checking indeed that 5(3;^'^°™'™'^^) (m 

^geom,conj (n-) • 

The two critical points of which 
correspond to the critical points of V 
(as ti — > 0) are indicated in Figures |^ 
and ^. As mentioned at the end of Sec- 



,m 



V ^ X 

V ^ X 

y X X A 

X 

V ^ X 
X 

X y 

5t Y^x^ 
■^Kx X 

■X X 
.X 



-15 -10 -5 



tion 3.4, at any fixed h there exist 
many other critical points of Afi that 
can be seen between consecutive pairs 
of poles and zeroes in Figures |^ and 
|8|; however, these other critical points 
become trapped in half-line singulari- 
ties as h ^ 0, and their saddle point 
approximations are not well-defined. 
We now calculate the perturbative 



10 



15 



Re(p) 



n " O 
3 O On 

^ol^ -151 



Figure 9: Poles, zeroes, and critical points of e^'y^-v^'u) 
for M = and n = |e*^/^ 



invariants 5'i^°°™''(n) and Sn"^^''\u) by doing a full saddle point approximation of the in- 
tegral ( [4.46| ) on (dummy) contours passing through the two critical points. We begin by 
^^Here s(a;,m) is the square root of a rational function rather than a rational function itself due to our 



redefined variables. Using (4.4E) directly, we would have gotten a pure rational expression. 
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formally expanding T{h,p,u) = logAfi{p,u) as a series in both h and p around some fixed 
point pq: 

oo oo 

(4.51) 



T{h,po +p,u) = J2 J2 '^j,kipo,u)^h^ 

j=0 fc=-l 

Om' potential V{p,u) is identified with To,-i(p, u) + liriu, and critical points of V are 
defined by Ti__i(p, u) = -^V{p,u) = 0. Let us also define 

6(p,^/) := -2T2-i{p,u) = —^V{p,u) . 
Then at a critical point po = p^"\u), the integral ( 4.46| ) becomes^^ 



(4.52) 



Z'^°'\M;h,u) = 



„«+-iy<=)(u) 



c?pe 2n P exp 



Co 



00 00 



3=3 



3=0 k=i 



where V^'^\u) = V {p^'^\u) , u) , 6(")(u) = 6(p(°)(u), ^x), and Tj."^^(^i) = T^- fc(p(")(n), u) 
implicitly functions of u alone. 

We can expand the exponential in ( [4.53 ), integrate each term using 



(4.53) 
are 



/ 



dpe~ 



b „2 



(n-l)!!(f)"/^y^ neven 
n odd 



(4.54) 



and re-exponentiate the answer to get a final result. The integrals in ( [4.54 ) are accurate 
up to corrections of order 0(6"'=°°^*/'^)^ which depend on a specific choice of contour and 
are ignored. Following this process, we obtain 



1 



V2^ V 6(°) 
exp 



27r/i ^+j.v('^)(u) sirh+s^^"h^+.. 



V(-)(n) - - log + u + si''h + s^^h^ + ... 



(4.55) 
(4.56) 



where the coefficients Sn can be straightforwardly computed in terms of 6 and the T's. For 



(a) 



15 "r(") I 3 "r(") I t(°^) cr^A q(") 



example, 

other terms). In addition, we clearly have 

5j")(u) = y(°)(n), 

= , 



1 

o log 



2 • 



^(t(^,)^+ (sixteen 

(4.57) 
(4.58) 

(4.59) 



To actually evaluate the coefficients Tij{p, u), we refer back to the expansion ( [3. 561 ) of the 
quantum dilogarithm in Section We find 



(4.60) 



T„fc(p,n) = ^'^I'l^l^^f [Lii_,_fc(-e^-") - (-l)^Lii„,„,(-e-^-'^) 



^This expression assumes that Tj,k ~ when k is even, a fact that shall be explained momentarily. 
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when J > 2 or A; > 0, and that all '^j,2k vanish. 

Therefore: to calculate the expansion coefficients sli^\u) around a given critical point, 
we substitute p^"\u) from ( 4.50 ) into equations (|4.48| ) and ( |4.60 ) to obtain b^°'\ and 

T^); then we substitute these functions into expressions for the Sn and simplify. At the 
geometric critical point p^^°"^\u), we obtain 



(4.61) 

2 

6feeom)(^) ^ !!ILA(m), (4.62) 



and it is easy to check with a little algebra that all the expansion coefficients si^^"^\u) 
reproduce exactly^^ what we found in the previous subsection by quantizing the moduli 
space of flat connections. (This has been verified to eight-loop order.) Similarly, at the 
conjugate critical point p(™'^j)(u), we have 



2 

y(=°^j)(n) = -y(s«°°^)(-u) , 6(™''j)(u) = -^A(m) , (4.63) 

and more generally Sn°^^^ = (— l)"'~"'^5'i^°°™\ It is not hard to actually prove this re- 
lation between the geomtric and conjugate critical points to all orders by inspecting the 
symmetries of e^^^'^'^\ Thus, we find complete agreement with the results of Section 4.2.1 



Acknowledgments 

We would like to thank D. Auroux, N. Dunfield, S. Garoufalidis, K. Hikami, T. Mrowka, 
W. Neumann, E. Witten, and C. Zickert for useful discussions and correspondence. Re- 
search of SG is supported in part by NSF Grant PHY-0757647 and in part by the Alfred P. 
Sloan Foundation. JL acknowledges support from a Marie Curie Intra-European Fellow- 
ship. TD acknowledges support from a National Defense Science and Engineering Graduate 
Fellowship. Opinions and conclusions expressed here are those of the authors and do not 
necessarily reflect the views of funding agencies. 



References 

[1] J.E. Andersen and S.K. Hansen, Asymptotics of the quantum invariants for surgeries on the 
figure 8 knot, J. Knot Theory and Its Ramifications 15 (4), 479-548 (2006). 

[2] M. Atiyah, The Geometry and Physics of Knots, Cambridge Univ. Press, 1990. 

[3] S. Axelrod and I. M. Singer, Chern-Simons perturbation theory, Proceedings of the XXth 
International Conference on Differential Geometric Methods in Theoretical Physics, Vol. 1, 2 
(New York, 1991) (River Edge, NJ), World Sci. Publishing, 1992, pp. 3-45. 



^'^There appears to be a small "correction" of ilog(— 1) = log(±j) in g^scom)^ comparing (4.5S:) with 
the value in Table ^. This merely multiplies the partition function by i and can be attributed to the 
orientation of the stationary-phase contour passing through the geometric critical point. We also allow the 
usual modulo 2niu ambiguity in matching 5*0. 



-57- 



[4] D. Bar-Natan, Pcrturhativc Aspects of the Chcrn-Simons Topological Quantum Field 
Theory, Ph.D. thesis, Princeton Univoristy, Juno 1991. 

[5] D. Bar-Natan and E. Witten, Perturbative expansion of Chern-Simons theory with 
noncompact gauge group, Commun. Math. Phys. 141, 423 (1991). 

[6] F. Bayen, M. Flato, C. Pronsdal, A. Lichnerowicz and D. Sternheimer, Deformation Theory 
And Quantization. 1. Deformations Of Symplectic Structures, Annals Phys. Ill, 61 (1978). 

[7] S. Bloch, H. Esnault and D. Kreimer, On motives associated to graph polynomials, 
Commun. Math. Phys. 267, 181-225 (2006), math/0510011. 

[8] D. Cooper, M. Culler, H. Gillet, D. Long and P. Shalen, Plane curves associated to character 
varieties of 3-manifolds, Invent. Math. 118, 47-84 (1994). 

[9] D. Cooper and D. Long, Remarks on the A-polynomial of a knot, J. Knot Theory and Its 
Ramifications 5 (5), 609-628 (1996). 

[10] O. Costin and S. Garoufalidis, Resurgence of 1-dimensional sums of Sum-Product type, in 
preparation. 

[11] R. Dijkgraaf and H. Fuji, The VbJume Conjecture and Topological Strings, arXiv:0903.2084 
[hep-th]. 

[12] T. Dimofte and J. Lenells, Observations on the State Integral Model for SL{2, C) 
Chern-Simons Theory, in preparation. 

[13] L. D. Faddeev, Current-Like Variables in Massive and Massless Integrable Models, in 
Quantum groups and their applications in physics (Varenna 1994), pages 117-135, Proc. 
Internat. School Phys. Enrico Fermi, 127, lOS, Amsterdam, 1996, hep-th/9408041. 

[14] L. Faddeev, Modular Double of Quantum Group, (1999), math/9912078. 

[15] L. D. Faddeev, R. M. Kashaev and A. Y. Volkov, Strongly coupled quantum discrete 

Liouville theory. I: Algebraic approach and duality, Commun. Math. Phys. 219, 199-219 
(2001), hep-th/0006156. 

[16] D. S. Freed, Remarks on Chern-Simons Theory, arXiv:0808.2507 [math. AT]. 

[17] S. Garoufalidis, On the characteristic and deformation varieties of a knot, in Geometry and 
Topology Monographs, number 7, pages 291-309, 2004. 

[18] S. Garoufalidis and J. Geronimo, Asymptotics of q-difference equations, 2006. 

[19] S. Garoufalidis and J. Geronimo, A Riemann-Hilbert approach to the asymptotics of 
1-dimensional sums of Sum-Product type, in preparation. 

[20] A.B. Goncharov, The pentagon relation for the quantum dilogarithm and quantized Mo,5, in 
Geometry and dynamics of groups and spaces, pages 415-428, Progr. Math. 265, Birkhauser, 
Basel, 2008, arXiv:0706.4054 [math.QA]. 

[21] D. Griinberg and P. Moree, Sequences of enumerative geometry: congruences and 
asymptotics. Experim. Math., to appear. 

[22] S. Gukov, Three-Dimensional Quantum Gravity, Chern-Simons Theory, and the 
A-Polynomial, Commun. Math. Phys. 255, 577-627 (2005), hep-th/0306165. 

[23] S. Gukov and H. Murakami, SL(2,C) Chern-Simons theory and the asymptotic behavior of 
the colored Jones polynomial, (2006), math/0608324. 



-58- 



[24] S. Gukov and E. Witten, Branes and Quantization, arXiv:0809.0305 [hep-th]. 

[25] O. Heinonen, (cd.) Composite fermions: A unified view of the quantum Hall regime, 
Singapore: World Scientific, 1998. 

[26] K. Hikami, Hyperbolic Structure Arising from a Knot Invariant, Int. J. Mod. Phys. A 16, 
3309-3333 (2001), math-ph/0105039. 

[27] K. Hikami, Generalized Volume Conjecture and the A-Polynomials - the Neumann- Zagier 
Potential Function as a Classical Limit of Quantum Invariant, 2006. 

[28] N. Hitchin, The Self-Duality Equations On A Riemann Surface, Proc. London Math. Soc.(3) 
55, 59-126 (1987). 

[29] R. M. Kashaev, The hyperbolic volume of knots from the quantum dilogarithm, Lett. Math. 
Phys. 39 269 275 (1997), q-alg/9601025. 

[30] R. M. Kashaev, On the spectrum of Dehn twists in quantum Teichmuller theory, in Physics 
and combinatorics (Nagoya, 2000), pages 63-81, World Sci. Publ., River Edge, NJ, 2001, 
math/0008148. 

[31] M. Kontsevich and D. Zagier, Periods, in Mathematics unlimited — 2001 and beyond, pages 
771-808, Springer, BerHn, 2001. 

[32] M. Kontsevich, Deformation Quantization Of Algebraic Varieties, Lett. Math. Phys. 56 
271-294 (2006), math/0106006. 

[33] J. M. F. Labastida, Chcrn-Simons gauge theory: Ten years after, (1999), hep-th/9905057. 

[34] J. Lewis, D. Zagier, Period functions for Maass wave forms. I, Ann. of Math. 153, 191-258 
(2001). 

[35] M. Marino, Chern-Simons theory and topological strings, Rev. Mod. Phys. 77 (2005) 675, 
hep-th/0406005. 

[36] J. Milnor, Hyperbolic geometry: the first 150 years. Bull. Amer. Math. Soc. 6, 9-24 (1982). 

[37] CD. Mostow, Strong rigidity of locally symmetric spaces, Ann. Math. Studies 78 (1973). 

[38] H. Murakami, A version of the volume conjecture, (2006), math/0603217. 

[39] G. Murthy and R. Shankar, Hamiltonian theories of the fractional quantum Hall effect. Rev. 
Mod. Phys. 75 (2003) 1101. 

[40] W. D. Neumann, Extended Bloch group and the Cheeger-Chern-Simons class, Geom. Topol. 
8, 413 474 (2004). 

[41] W. D. Neumann and J. Yang, Bloch invariants of hyperbolic 3-manifolds, Duke Math. J. 

96, 29 59 (1999). 

[42] W. D. Neumann and D. Zagier, Vbiumes of hyperbolic three-manifolds, Topology 24, 
397-332 (1985). 

[43] C. Petronio and J.R. Weeks, Partially fiat ideal triangulations of cusped hyperbolic 
3-manifolds, Osaka J. Math. 37, 453-466 (2000). 

[44] G. Prasad, Strong rigidity of Q -rank 1 lattices, Invent. Math. 21, 255-286 (1973). 

[45] N. Seiberg and E. Witten, Electric - magnetic duality, monopole condensation, and 
confinement in N=2 supersymmetric Yang-Mills theory, Nucl. Phys. B 426 (1994) 19, 
[Erratum-ibid. B 430 (1994) 485], hep-th/9407087. 



-59- 



[46] W. P. Thurston, The geometry and topology of three-manifolds, Lecture notes at Princeton 

University, Princeton, 1980. 

[47] W. P. Thurston, Three Dimensional Manifolds, Kleinian Groups, and Hyperbolic Geometry, 
Bull. AMS 6 (3), 357-381 (1982). 

[48] W. P. Thurston, Hyperbolic structures on 3-manifolds. I. Deformations of acylindrical 
manifolds, Ann. Math. 124 (2), 203-246 (1986). 

[49] D. Thurston, Hyperbolic Volume and the Jones Polynomial, in Invariants de Noeuds et de 
Varietes de Dimension 3 (Greboble 1999), Institut Fourier, 1999. 

[50] V.G. Turaev and O.Yu. Viro, State-sum invariants of 3-manifolds and quantum 6j-symbols, 
Topology 31 865 (1992). 

[51] A. Y. Volkov, Nonconunutative Hyper geometry, Commun. Math. Phys. 258 (2), 257-273 
(2005), math/0312084. 

[52] E. Witten, Quantum field theory and the Jones polynomial, Commun. Math. Phys. 121, 
351-399 (1989). 

[53] E. Witten, Quantization Of Chern-Simons Gauge Theory With Complex Gauge Group, 
Commun. Math. Phys. 137, 29 (1991). 

[54] N. Woodhouse, Geometric Quantization, New York: Oxford University Press, 1992. 

[55] T. Yoshida, The rj-invariant of hyperbolic 3-manifolds, Invent. Math. 81, 473-514 (1985). 

[56] D. Zagier, Vassiliev invariatnts and a strange identity related to the Dedekind eta-function, 
Topology 40, 945-960 (2001). 

[57] D. Zagier, The dilogarithm function, In Frontiers in Number Theory, Physics and 
Geometry II, P. Carticr, B. Julia, P. Moussa, P. Vanhove (eds.). Springer- Verlag, 
Berlin-Heidelberg-New York, 3 65 (2006). 

[58] D. Zagier, Algebraic and asymptotic properties of quantum invariants of knots, in 
preparation. 

[59] D. Zagier, Quantum modular forms. Preprint (2009), 16 pages. 

[60] C. Zickert, The Chern-Simons invariant of a representation, arXiv:0710.2049 [math.GT]. 



-60- 



